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TIME VARYING SOLAR CYCLE PROTONS 
PROGRAM MANUAL 


INTRODUCTION 


The purpose of this manual is to provide a detailed description of the pro- 
grams employed in calculating the variation of the proton population in the earth's 
radiation belt due to solar cycle atmospheric and source changes. This report 
makes no attempt to make detail explanations of the theory involved in the cal- 
culations. An adequate list of references is provided for those interested in the 
explanation of the methods and the results to be expected. The entire process 
involves several interdependent steps. For example, the diurnal averaged at- 
mosphere and the rings output from the B-L search routine are input to the 
longitudinal averaging processor. The output from this program is then used as 
input to the "bounce” average program and so forth. Each program description 
includes the equations used, a flow chart and Fortran listing of the program, 
input and output specifications, card descriptions, sample input and output data 
and running time. The code is written in Fortran n language and assembled under 
GSFC-NASA Theoretical Division monitor system on the IBM-7094 computer. 
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I. DIURNAL AVERAGED ATMOSPHERE 


Contained here are tables of number density vs. altitude and solar flux num- 
ber. Separate tables of He, 0, 0 2 , N 2 and H are contained in Tables 1-5. The 
models are those generated by Harris and Priester . ( 2 ) Each model refers to 
a given solar radiation flux in units of 10- 22 watts/m 2 /cycle/sec. The link be- 
tween solar flux, S, and time is given by Figure 1 (reproduced from reference 2 
with the permission of Harris and Priester). This data is included for continuity 
and is used as input for the longitudinal averaging processor. The Harris and 
Priester atmosphere and solar flux vs. time curve is the data used and is subject 
to change as soon as better data becomes available. Such changes can be per- 
formed without effecting the logic of the following programs. 
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Table 1 

Diurnal Averaged Number Densities of He as a 
Function of Altitude for Five Solar Flux Numbers. 


s 

h(krru5 s v 

250 

200 

150 

100 

70 

120 

2.500E07 

2.500E07 

2.500E07 

2.500E07 

2.500E07 

200 

4.413E06 

4.872E06 

5.484E06 

6.319E06 

6.982E06 

300 

2.642E06 

2.843E06 

3.050E06 

3.205E06 

3.212E06 

400 

1.945E06 

2.007E06 

2.012E06 

1.885E06 

1.689E06 

500 

1.499E06 

1.477E06 

1.390E06 

1.154E06 

9.272E05 

600 

1.176E06 

1.107E06 

9.646E05 

7.236E05 

5.231E05 

700 

9.327E05 

8.387E05 

6.842E05 

4.607E05 

3.021E05 

800 

7.461E05 

6.422E05 

4.916E05 

3.004E05 

1.783E05 

900 

6.016E05 

4.964E05 

3.574E05 

1.984E05 

1.074E05 

1000 

4.886E05 

3.871E05 

2.628E05 

1.331E05 

6.593E04 

1100 

3.995E05 

3.044E05 

1.953E05 

9.050E04 

4.119E04 

1200 

3.287E05 

2.411E05 

1.466E05 

6.238E04 

2.616E04 

1300 

2.721E05 

1.925E05 

1.110E05 

4.353E04 

1.687E04 

1400 

2.265E05 

1.547E05 

8.482E04 

3.074E04 

1.104E04 

1500 

1.895E05 

1.251E05 

6.537E04 

2.195E04 

7.328E03 

1600 

1.595E05 

1.019E05 

5.078E04 

1.584E04 

4.925E03 

1700 

1.348E05 

8.346E04 

3.974E04 

1.155E04 

3.351E03 

1800 

1.145E05 

6.875E04 

3.133E04 

8.497E03 

2.306E03 

1900 

9.774E04 

5.696E04 

2.487E04 

6.308E03 

1.605E03 

2000 

8.378E04 

4.743E04 

1.988E04 

4.724E03 

1.128E04 







Table 2 

Diurnal Averaged Number Densities of O as a 
Function of Altitude for Five Solar Flux Numbers. 


\ s 

h(k^u^'v 

250 

200 

150 

100 

70 

120 

7.600E10 

7.600E10 

7.600E10 

7.600E10 

7.600E10 

200 

3.600E09 

3.457E09 

3.209E09 

2.795E09 

2.416E09 

300 

8.870E08 

7.134E08 

5.124E08 

2.809E08 

1.564E08 

400 

3.054E08 

2.054E08 

1.112E08 

4.025E07 

1.471E07 

500 

1.168E08 

6.616E07 

2.788E07 

6.771E06 

1.675E06 

600 

4.749E07 

2.287E07 

7.708E06 

1.273E06 

2.183E05 

700 

2.024E07 

8.364E06 

2.232E06 

2.611E05 

3.153E04 

800 

8.983E06 

3.207E06 

6.918E05 

5.747E04 

4.946E03 

900 

4.130E06 

1.282E06 

2.252E05 

1.342E04 

8.320E02 

1000 

1.960E06 

5.312E05 

7.645E04 

3.301E03 

1.488E02 

1100 

9.567E05 

2.275E05 

2.696E04 

8.502E02 

2.810E01 

1200 

4.791E05 

1.003E05 

9.834E03 

2.284E02 

5.381E00 

1300 

2.456E05 

4.538E04 

3.701E03 

6.379E01 

1.162E00 

1400 

1.287E05 

2.106E04 

1.434E03 

1.848E01 

2.527E-1 

1500 

6.878E04 

9.997E03 

5.706E02 

5.541E00 

5.729E-2 

1600 

3.746E04 

4.849E03 

2.330E02 

1.717E00 

1.352E-2 

1700 

2.077E04 

2.399E03 

9.744E01 

5.486E-1 

3.312E-3 

1800 

1.170E04 

1.210E03 

4.170E01 

1.806E-1 

8.412E-4 

1900 

6.700E03 

6.216E02 

1.824E01 

6.079E-2 

2.212E-4 

2000 

3.893E03 

3.247E02 

8.149E00 

2.128E-2 

6.01 0E-5 
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Table 3 

Diurnal Averaged Number Densities of 0 2 as a 
Function of Altitude for Five Solar Flux Numbers. 



1.200E11 

9.900E08 

7.683E07 

1.020E07 

1.217E06 

3.142E05 

6.428E04 

1.364E04 

3.242E03 

7.981E02 

2.056E02 

5.526E01 

1.544E01 

4.474E00 

1.342E00 

4.157E-1 

1.328E-1 

4.370E-2 

1.478E-2 

5.142E-3 


1.200E11 

7.910E08 

4.269E07 

4.069E06 

4.886E05 

6.758E04 

1.033E03 

1.714E03 

3.031E02 

5.682E01 

1.122E01 

2.322E00 

5.028E-1 

1.136E-1 

2.778E-2 

6.520E-3 

1.653E-3 

4.342E-4 

1.177E-4 

2.906E-5 


1.200E11 
5.699E08 
1.791E07 
1.048E06 
7.983E04 
7.106E03 
7.062E02 
7.688E01 
8.962E00 
1.122E00 
1.493E-1 
2.106E-2 
3.139E-3 
4.927E-4 
8.124E-5 
1.405E-5 
2.541E-6 
4.800E-7 
9.452E-8 
1.937E-8 











Table 4 


Diurnal Averaged Number Densities of N 2 as a 
Function of Altitude for Five Solar Flux Numbers. 


s 

h(krm)\. 

250 

200 

150 

100 

70 

120 

5.800E11 

5.800E11 

5.800E11 

5.800E11 

5.800E11 

200 

7.393E09 

6.180E09 

4.743E09 

3.136E09 

2.124E09 

300 

7.630E08 

4.639E08 

2.210E08 

6.739E07 

2.113E07 

400 

1.278E08 

5.777E07 

2.151E07 

2.682E06 

4.225E05 

500 

2.562E07 

8.798E06 

1.810E06 

1.407E05 

1.156E04 

600 

5.730E06 

1.521E06 

3.239E05 

8.800E03 

3.844E02 

700 

1.403E06 

2.884E05 

2.763E04 

6.190E02 

1.471E01 

800 

3.664E05 

5.883E04 

3.899E03 

4.802E01 

6.320E-1 

900 

1.012E05 

1.276E04 

5.894E02 

4.057E00 

3.006E-2 

1000 

2.932E04 

2.918E03 

9.470E01 

3.702E-1 

1.570E-3 

1100 

8.862E03 

6.994E02 

1.609E01 

3.627E-2 

8.936E-5 

1200 

2.783E03 

1.725E02 

2.878E00 

3.798E-3 

1.016E-5 

1300 

9.053E02 

4.558E01 

5.404E-1 

4.237E-4 

3.674E-7 

1400 

3.041 E02 

1.233E01 

1.063E-1 

5.017E-5 

2.633E-8 

1500 

1.075E02 

3.456E00 

2.183E-2 

6.287E-6 

2.023E-9 

1600 

3.768E01 

1.002E00 

4.677E-3 

8.321E-7 

1.661E-10 

1700 

1.380E01 

3.001E-1 

1.042E-3 

1.160E-7 

1.454E-11 

1800 

5.211E00 

9.268E-2 

2.313E-4 

1.699E-8 

1.353E-12 

1900 

2.O05EOO 

2.949E-2 

5.746E-5 

2.612E-9 

1.334E-13 

2000 

7.927E-1 

9.654E-3 

1.431E-5 

4.619E-10 

1.393E-14 
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Table 5 


Diurnal Averaged Number Densities of H as a 
Function of Altitude for Five Solar Flux Numbers. 


s 

h(kmO N v 

250 

200 

150 

100 

70 

120 

4.356E04 

4.356E04 

4.356E04 

4.356E04 

4.356E04 

200 

1.071E04 

1.224E04 

1.447E04 

1.790E04 

2.104E04 

300 

8.035E03 

9.323E03 

1.114E04 

1.380E04 

1.611E04 

400 

7.205E03 

8.328E03 

9.837E03 

1.189E04 

1.352E04 

500 

6.690E03 

7.660E03 

8.898E03 

1.046E04 

1.157E04 

600 

6.272E03 

7.102E03 

8.107E03 

9.263E03 

9.974E03 

700 

5.904E03 

6.609E03 

7.415E03 

8.214E03 

8.638E03 

800 

5.573E03 

6.168E03 

6.803E03 

7.360E03 

7.518E03 

900 

5.272E03 

5.769E03 

6.260E03 

6.597E03 

6.571E03 

1000 

4.996E03 

5.408E03 

5.774E03 

5.933E03 

5.768E03 

1100 

4.742E03 

5.079E03 

5.339E03 

5.352E03 

5.083E03 

1200 

4.502E03 

4.778E03 

4.947E03 

4.843E03 

4.495E03 

1300 

4.291E03 

4.503E03 

4.594E03 

4.395E03 

3.989E03 

1400 

4.090E03 

4.250E03 

4.275E03 

3.998E03 

3.552E03 

1500 

3.903E03 

4.018E03 

3.986E03 

3.647E03 

3.173E03 

1600 

3.730E03 

3.804E03 

3.723E03 

3.335E03 

2.843E03 

1700 

3.568E03 

3.606E03 

3.484E03 

3.057E03 

2.555E03 

1800 

3.417E03 

3.424E03 

3.266E03 

2.809E03 

2.303E03 

1900 

3.276E03 

3.255E03 

3.066E03 

2.587E03 

2.081E03 

2000 

3.144E03 

3.098E03 

2.884E03 

2.387E03 

1.886E03 
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H. B-L SEARCH 


A. Introduction 

This program produces contours of constant B and L as a function of longi- 
tude, latitude and altitude in both the northern and southern hemispheres. The 
desired initial values of B and L are read into the program along with an 
approximate corresponding latitude, K , which can easily be obtained by the use 
of Figure 2. The dipole equation r Q = L 0 cos 2 A. relating the initial L and 
the geocentric distance r Q is used with the radius of the earth, r e , and the equa- 
tion hj = r e (r Q - 1) to provide an approximation of the altitude. This is fed, 
together with a longitude of 18C degrees and latitude k, into subroutine INVAR 
which calculates B and L for a given longitude, latitude and altitude. This sub- 
routine makes use of the transformation developed by Mcllwain ( 3 using the 48 
spherical harmonic coefficients of Jenson and Cain (4) . INVAR numerically 
integrates the longitudinal invariant I using a series expansion for the magnetic 
field. Then L is calculated as a function of both B and I by using a dipole rep- 
resentation of the earth. The B and L obtained in this manner are returned to 
the main program. Here the accuracy of the initial approximation is checked. 

If the computed L is found to be within an accuracy of 10 -4 of the initial L the 
program will enter into a search routine with linear interpolation in latitude 
and altitude in order to arrive at a correct B. The search parameters Ah (in- 
crement in altitude) and A A. (increment in latitude) are prefixed and must remain 
small in order that interpolation may hold. Once B is found, it is checked, to- 
gether with the value of L, to insure an accuracy of 10~ 3 in comparison with the 
initial values. If the accuracy is sufficient the subroutine RING will be called 
for the northern hemisphere. This subroutine takes a given latitude and altitude 
and computes the B and L contour map for longitudes of 10 degree increment 
for the full 360 degrees. The program then computes the contour map for the 
southern hemisphere at the same B and L, increments B by .01 gauss and returns 
to the northern hemisphere. It will continue in this manner until the altitude 
drops below 100 kilometers, at which point the next initial B, L and k will be 
read until input data is exhausted. The entire process prepares input for the 
longitudinal averaging processor. 
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B. Mnemonics 




Quantity 

Description 

Units 

CONVR 

conversion factor — degrees to radians 

— 

RE 

radius of the earth 

km. 

DELH 

delta altitude 

km. 

ALONG 

geocentric longitude 

degrees 

DLAT 

delta latitude 

tt 

ALAT 

geocentric latitude 

TT 

B 

magnetic induction 

gauss 

EL 

magnetic field line 

earth radii 

R 

geocentric distance 

TT 

H 

altitude 

km. 

SAVB 

temporary storage of B 

gauss 

SAVH 

temporary storage of H 

km. 

SAVLAT 

temporary storage of ALAT 

degrees 

LCNT 

an indicator used to indicate when 
the altitude has dropped below 100 km. 

- - 
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E. Restrictions 


ALAT must be remotely close to the correct latitude corresponding to the 
initial B in order to save machine time and to insure the accuracy of the answers. 
ALAT must always be positive and greater than the geomagnetic equator at geo- 
centric longitude of 180°. This means that ALAT must be greater than 4°. 

F . Input 

Cards containing the initial values of B, L and \ are all the data necessary 
for the execution of this program. This data is entered on logical tape number 
two. Each card represents a single case. 


1. Card Description 


Columns 

Mode 

Quantity 

Units 

Description 

1-12 

F 

BO 

gauss 

initial magnetic induction 

13-24 

F 

ELO 

earth radii 

initial magnetic field line 

25-32 

F 

ALAT 

degrees 

initial latitude 
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GSPC FORM 541-1 (July - 40) 



G. Output 


Output from this program appears on logical tapes 3 and 5. Both tapes 
contain the same information with the exception that tape 3 also contains the 
intermediate results of the search routine. Tape 5 is the output to be used in 
the longitudinal averaging processor. This tape is punched onto cards which 
can be combined with similar output from other runs. All this data can then be 
used as input to the longitudinal averaging processor, but it is important that a 
blank card follow the last of the data cards. This arrangement is discussed in 
the next section. 
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1. Tape 3 Sample 


B0= 0.11C0Q L0= 1.75000 INITIAL LAT=26. 00000 

INTERMEDIATE RESULTS — NORTHERN HEMISPHERE 


LONG 

LAT 

ALT 

8 

L 

180.C0CC0 

26.00000 

3533.22839 

0.10159 

1 . 75000 

180.00000 

28.00000 

3257.93388 

0.11298 

1.75000 

180.00000 

27.87628 

3330.02185 

0.10985 

1.74955 

0.10985 

L= 1.78955 






FINAL RESULTS 



-180.00 

27.88 

3330.02 

0.10985 

1.74955 

-180.00 

27.88 

3330.02 

0.10965 

1 . 74955 

-1 70.00 

25.79 

3309.07 

0. 10987 

1.75027 

-160.00 

23.86 

3298.23 

0.10987 

1.75043 

-150.00 

21. 82 

3293.89 

0.10987 

1.75049 

-180.00 

19.78 

3290.96 

0.10986 

1.74966 

-130.00 

17.71 

3287.81 

0. 1C986 

1.74966 

-120.00 

15.73 

3279.11 

0. 10986 

1.74966 

-110.00 

13.83 

3263.27 

0.10987 

1.75054 

-ICO. 00 

11.97 

3238.87 

0. 10987 

1.75049 

-90.00 

10.27 

3190.87 

0.10987 

1.75052 

-80.00 

8.93 

3132.92 

0.10986 

1 .75014 

-70.00 

8.27 

3067.73 

0.10986 

1.74978 

-60.00 

8.66 

3009.01 

0.10985 

1.74947 

-50.00 

10. 36 

2972.82 

0. 10985 

1.74926 

-80.00 

13. 12 

2968.85 

0.10985 

1.74926 

-30.00 

16. 35 

2979.89 

0.10985 

1.74933 

-20.00 

19.88 

3011.16 

0.10985 

1.74939 

-10.00 

22. 19 

3052.95 

0.10985 

1.74927 

-0. 

28.38 

3101.38 

0.10985 

1.74913 

10.00 

25.90 

3153.89 

0.10965 

1.74914 

20.00 

26.92 

3210.03 

0.10985 

1.74922 

30.00 

27.89 

3269.30 

0.10985 

1.74940 

80.00 

27.80 

3331.19 

0.10986 

1.74874 

50.00 

28.10 

3397.68 

0.10987 

1.74857 

60.00 

28.53 

3866.52 

0.10985 

1.74941 

70.00 

29.03 

3528.77 

0.10985 

1.74939 

80.00 

29.55 

3579.90 

0.10985 

1.74936 

90.00 

29.99 

3618.76 

0.10985 

1.74941 

100.00 

30.2 7 

3629.99 

0.10985 

1.74857 

110.00 

30.88 

3626.67 

0.10985 

1.74880 

120.00 

30.6 3 

3605.07 

0.10986 

1.74932 

130.00 

30.6 7 

3566.83 

0. 10987 

1.74938 

180.00 

30.68 

3517.70 

0. 10987 

1.74968 

150.00 

30.82 

3868.00 

0.10988 

1.75053 

160.00 

29.83 

3810.37 

0.10986 

1.74985 

170.00 

28.87 

3368.67 

0. 10986 

1.75002 
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INTERMEDIATE RESULTS — SOUTHERN HEMISPHERE 


LCM3 

LAT 

ALT 

L' 

L 

ISO. CCCGO 

— 15.000 Cu 

3695.08514 

0. 10030 

1 . 74996 

180.00000 

-17.00000 

3433.09256 

0.11171 

1. 74995 

180.00000 

-16.70013 

3472.37378 

0.10990 

1.74989 

C. 10990 

L = 1.7498 

9 





FINAL RESULTS 



-160.00 

-16.70 

3472.37 

0.10990 

1.74989 

-180.00 

-16.70 

3472.37 

0.10990 

1.74989 

-170.00 

-18.67 

3441.67 

0.10990 

1.75006 

-160.00 

-20.70 

3407.85 

0.10990 

1.75001 

-150.00 

- 22.78 

3370.58 

0.10990 

1.75002 

1 

-C" 

o 

• 

o 

O 

-24.86 

3327.75 

0.10990 

1.75005 

-130. CO 

-26.87 

3277.04 

0.10988 

1.74890 

-120.00 

-28.92 

3219.97 

0.10987 

1 . 74896 

-110.00 

-30.96 

3153.88 

0.10987 

1 . 74892 

-100.00 

-32.95 

3077.77 

0.10987 

1.74908 

-90.00 

-34.76 

2990.94 

0.10987 

1.74901 

-80.00 

-36.23 

2895.15 

0. 10966 

1.74917 

-70.00 

-37.17 

2794.22 

0.10989 

1.74943 

-60.00 

-37.39 

2694.25 

0.10990 

1.74975 

-50.00 

-36. 73 

2603.75 

0.10991 

1. 75009 

-40.00 

-35.09 

2533.92 

0.10991 

1. 75021 

-30.00 

-32.47 

2497.69 

0.10991 

1.75013 

-20.00 

-29.06 

2503.03 

0.10987 

1.74902 

-10.00 

-25.47 

2545.84 

0.10990 

1.74915 

-0. 

- 22.21 

2615.13 

0.10990 

1.74996 

10.00 

-19.60 

2695.72 

0.10990 

1.75003 

20.00 

-17.80 

2780.87 

0.10990 

1.75013 

30.00 

-16.70 

2870.66 

0.10990 

1 . 75012 

40.00 

-15.98 

2969.88 

0. 10990 

1.75012 

50.00 

-15.29 

3079.53 

0. 1 1000 

1. 74901 

60.00 

-14.53 

3205.88 

0.10990 

1. 75006 

70.00 

-13.64 

3328.53 

0.10990 

1.75010 

80.00 

-12.77 

3437.57 

0.10990 

1.75010 

90.00 

-12.08 

3524.42 

0. 10990 

1.75011 

100.00 

-11.65 

3585.00 

0.10990 

1. 75007 

110.00 

-11.42 

3618.98 

0.10990 

1.74908 

120.00 

-11.44 

3631.77 

0.10990 

1 . 74995 

130.00 

-11.60 

3624.71 

0.10991 

1.75053 

140.00 

-11.89 

3601.71 

0.10990 

1 . 74967 

150.00 

-12.51 

3570.65 

0. 10989 

1.74957 

160.00 

-13.51 

3536. 5 5 

0.10989 

1.74939 

170.00 

-14.92 

3503.21 

0.10988 

1.74915 



a. Card Description 


Columns 

Mode 

Quantity 

Units 

Description 

1-12 

F 

PRLNG 

degrees 

geocentric longitude 

13-24 

F 

PHI 

degrees 

geocentric latitude 

25-36 

F 

ALT 

km. 

altitude 

37-48 

F 

B 

gauss 

magnetic induction 

49-60 

F 

EL 

earth radii 

magnetic field line 
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GSFC FORM 541-1 (July - 60) 






H. Running Time 

This program takes about ten minutes for each 100 kilometers in altitude. 
The initial minimum altitude is required in order to estimate the number of 
kilometers for a particular run. Use Figure 2 and the initial values of B, L and 
\ to arrive at an estimate of the initial altitude h. Subtracting 800 kilometers 
from this value will give the approximate initial minimum altitude of the first 
ring. This is the altitude used in figuring the running time. 
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ID. LONGITUDINAL AVERAGING PROCESSOR 



A. Introduction 

Input to this program are the five diurnal average number density tables 
(Tables 1-5) and the B-L contour rings output from the B-L search routine. The 
tables are interpolated (extrapolated) in order to obtain a density value for every 
ten degrees of longitude in the B-L contours. The densities are then added to- 
gether for each of the five flux models, and the resulting sum is divided by 36 to 
arrive at the longitudinally averaged number density. This is done for the north- 
ern and southern hemispheres separately and then these values are added together 
and divided by two in order to obtain one number density for each B and L and 
each of the five flux models. These final values are used as input to the lambda 
punch program. 

B. Mnemonics 


Quantity 

Description 

Units 

H(J) 

J 01 altitude from the diurnal averaged 
atmosphere tables. Tables 1-5 

km. 

S1(J),. . S5(J) 

diurnal averaged atmosphere densities 
for the J* altitude and fluxes of 250, 
200, 150, 100 and 70 x 10" 22 watts/m 2 / 
cycle/sec from Tables 1-5 

atoms/cm 3 

Al,. . A5 

temporary storage of S1(J), . . .,S5(J) 

ft 

SS1(J),. . SS5(J) 

natural logarithms of S1(J), . . ,,S5(J) 

— 

AVHA 

sum of EL 

earth radii 

AVBN1 , . . AVBN5 

density summation for the five flux 
models 

atoms /cm 3 

LONG 

geocentric longitude 

degrees 

PHI 

geocentric latitude 

degrees 
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Quantity 

Description 

Units 

ALT 

altitude 

km. 

B 

magnetic induction 

gauss 

EL 

magnetic field line 

earth radii 

Y(I) 

temporary storage of SSK(J) at altitude 
I = J depending on whether K = 1, 2, 3, 
4 or 5 

_ _ 

HA 

temporary storage of altitude 

km. 

YA 

natural logarithm of density at altitude 
ALT 

— 

BARN 

density at altitude ALT 

atoms /cm 3 

BARNAV(K) 

temporary storage of BARN for flux 
model K 

IT 

JUNK 

temporary storage of K 

— 

AVERH1 

average longitudinal value of EL for 
only one hemisphere 

earth radii 

AVERN1, .... AVERN5 

averaged densities for the five flux 
models and one hemisphere 

atoms/cm 3 

FINI 

end of file trigger for tape 6 

— 

ELI 

EL for northern hemisphere 

earth radii 

EL2 

EL for southern hemisphere 

ft 

B1 

B for northern hemisphere 

gauss 

B2 

B for southern hemisphere 

n 
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Quantity 

Description 

Units 

EN1, . . EN5 

the first time this designation appears 
it represents the northern hemisphere 
densities for the five flux models. The 
second time it appears it represents 
the densities for the north and south 
averaged together. 

atoms/cm 

EN12, . . ., EN52 

averaged southern hemisphere densities 
for the five flux models 

t! 

ITT 

counter to keep track of which hemi- 
sphere is being considered 

— 

ITT! 

counter to keep track of which atmos- 
pheric constituent is being considered 

— 

XLA(ITTI) 

first six letters of the constituent 
name designated by ITU 

— 

XLB(ITTI) 

last two letters of the constituent 
name designated by ITU 

— 

TSA 

first six letters of constituent name 

— 

TSB 

last two letters of constituent name 

— — 
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C. Flow Chart 



Use log extrapolation to 
find log of density for 
Kth flux model at altitude 


Use log extrapolation to 
find log of density for 
Kth flux model at altitude 


yes 



YtJ) 1 

ALT. call it YA. 


ALT. call it YA. 

i 

r 


Find exp(YA) 
(density at 
altitude ALT) 


Use log interpolation with 
this altitude to find log of 
density for the Kth flu* 
model at altitude ALT. 
call it YA. 
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EaIlRML FORMULA NUMBERS WITH CCRKlSPCNC I NG INTERNAL FORMULA NUMBERS AND OCTAL LOCATIONS 
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E. Restrictions 


Input tables must have densities for altitudes starting at 120 kilometers and 
followed by each 100-kilometer level up to 2000 kilometers. If the tables are 
changed then the program must be modified. The tables must be input in the 
following order: 

(1) helium, 

(2) oxygen, 

(3) molecular oxygen, 

(4) nitrogen, 

(5) hydrogen. 

F . Input 

Input to this program is handled by two tapes. On tape 2 are entered the 
diurnal averaged atmosphere tables (Tables 1-5). These tables are punched on 
cards and placed one behind the other at the end of the Fortran deck. These 
tables must be input in the following order: helium, oxygen, molecular oxygen, 
nitrogen and hydrogen. This is necessary in order that the constitutuent names 
may be correctly punched on the output. It is also important that a blank card 
follow the last table in order to notify the program when it has reached the end 
of the file. 

On tape 5 are placed the B-L contours which were produced by the B-L 
search program. As with the tables, the last card must be blank in order to 
designate an end of file. 

1. Tape 2 


a. Input Card Description 


ATMOSPHERE TABLES 


Columns 

Mode 

Quantity 

Units 

Description 

1-12 

E 

S1(N) 

atoms/ cm 3 

density for flux model 1 , altitude N 

13-24 

E 

S2(N) 

tt 

Tt Tt 2 M 

25-36 

E 

S3(N) 

Tt 

n tt 2 n 

37-48 

E 

S4(N) 

Tt 

Tt TT 4 TT 

49-60 

E 

S5(N) 

IT 

" " 5, 
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GSFC FORM 541-1 (July - 60) 






2. Tape 5 


a. Input Card Description 




B-L 

CONTOURS 


Columns 

Mode 

Quantity 

Units 

Description 

1-12 

F 

LONG 

degrees 

geocentric longitude 

13-24 

F 

PHI 

degrees 

geocentric latitude 

25-36 

F 

ALT 

km 

altitude 

37-48 

F 

B 

gauss 

magnetic induction 

49-60 

F 

EL 

earth radii 

magnetic field line 
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GSFC FORM 541-1 (July - 60) 










G . Output 


Output for this program appears on logical tapes 3 and 10. Tape 3 will con- 
tain the averaged densities for each of the five flux models at a particular B and 
L for both northern and southern hemispheres and for each of the five constitu- 
ents. The constituents (He, 0, 0 2 , N 2 , H) follow each other in the same order as 
they exist in the tables which are input on tape 2. Tape 3 is printed with appro- 
priate headings. 

Tape 10 also contains densities for the five flux models at each B and L but 
here the northern and southern hemispheres are averaged together to give one 
final value for each B and L. Tape 10 contains no headings but each card is 
labeled with the appropriate constituent name. This tape is punched in order to 
be used as input to the lambda punch. 

A scratch tape must be set up on logical tape unit 6. Although no data ap- 
pears on this tape at the completion of the program it is used in intermediate 
steps to store data. 
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1. Tape 3 Sample 


ATMOSPHERIC CONSTITUENT — HELIUM 


AVERAGE L= 1.142 R=. 23891 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.8300b 06 0.7508E 06 0.6289E 06 0.4600E 06 0.3362E 06 


AVERAGE L= 1.142 8=. 23890 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.1346E 08 0.1210E 08 0.1064E 08 0.9134E 07 0.8207E C7 


AVERAGE L= 1.142 B=. 22741 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.7137E 06 0.6282E 06 0.5063E 06 0.3493E C6 0.2433E 06 


AVERAGE L= 1.142 8=. 22736 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

C.1428E 07 0.1402E 07 0.1357E 07 0.1293E 07 0.1250E 07 


AVERAGE L= 1.142 B=. 23465 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.7842E 06 0.7021E 06 0.5795E 06 0.4145E 06 0.2984b 06 


AVERAGE L= 1.14? B=. 23466 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.5205E 07 0.4932E 07 0.4608E 07 0.423SE 07 0.3986E 07 


AVERAGE L= 1.142 B=. 21968 

AVERAGE NUMBER DENSITIES FUR DIFF S 
NORTHERN HEMISPHERE 

0.6485E 06 0.5615E 06 0.4423E 06 0.2952E 06 0.1996E 06 


AVERAGE L= 1.142 B=. 21968 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.8665E 06 0.8125E 06 0.7325E 06 0.6263E 06 0.3477E Ob 


AVERAGE L - 1.170 B=. 23992 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.7300E 06 0.6416E 06 0.5L34E 06 0.3472E 06 0.2351E 06 


AVERAGE L= 1.170 B=. 23991 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.1571E 08 0.1399E 08 0.1217E 08 0.1032E 08 0.9197E 07 
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AVERAGE L = 1.170 B=. 23468 

AVERAGE NUMBER DENSITIES E OR DIFF S 
.NORTHERN HEMISPHERE 

0.6785E Oo 0.5881E 06 0.4614E 06 0.3026E 06 0.1990E 06 


AVERAGE L= 1.170 B=. 23467 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

C.5025E 07 0.4766E 07 0.4459E 07 0.4106E 07 0. 38 74E 07 


AVERAGE L= 1.170 B=. 22957 

AVERAGE NUMBER DENSITIES FUR DIFF S 
NORTHERN HEMISPHERE 

0.6319E 06 0. 5403t 06 0.4158E 06 0.2647E 06 0.1693E 06 


AVERAGE L = 1.170 B = . 22956 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.1975E 07 0.1940E 07 0.1885b 07 0.1810E 07 0.1758E C7 


AVERAGE L = 1.170 B = . 21976 
AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.5513E 06 0.4395E 06 0.3410E 06 0.2052E 06 0.1246E 06 


AVERAGE 1= 1.170 B=. 21977 

AVERAGE NUMBER DENSITIES FOR OIFF S 
SGUTHERN HEMISPHERE 

0.8559E 06 0.8062E 06 0.7343E 06 0.6406E 06 0.5740E C6 


AVERAGE L = 1.170 B=. 20975 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.4812E 06 0.3913E 06 0.2804E 06 0.1600E 06 0.9249E 05 


AVERAGE L= 1.170 B=. 20972 

AVERAGE NUMBER DENSITIES FUR DIFF S 
SOUTHERN HEMISPHERE 

0.6456E 06 0.5736E 06 0.4764E 06 0.3555E 06 0.2730E 06 


AVERAGE L= 1.L70 B=. 20474 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.4513E 06 0.3629E 06 0.2560E 06 0.1428E 06 0.8090E 05 


AVERAGE L = 1.170 B = . 20475 
AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.5692E 06 0.49UE 06 0.3896E 06 0.2688E 06 0.1507E 06 


AVERAGE L= 1.170 B=. 20141 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.4332E 06 0.3460E 06 0.2419E 06 0.1333E 06 0.7469E 05 
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AVERAGE L= 1.170 B=. 20134 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.5207E 06 0.4399E 06 0.3375E 06 0.2199E 06 0.1474E 06 


AVERAGE L= 1.188 B=. 23382 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.6179E 06 0.5248E 06 0.3992E 06 0.2485E 06 0.1550E 06 


AVERAGE L= 1.188 B=. 23383 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.37B2E 07 0.3632E 07 0.3447E 07 0.3229E 07 0.3084E 07 


AVERAGE E= 1.186 B=. 21177 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

C.4507E 06 0.3605E 06 0.2517E 06 0.1373E 06 0.7561E 05 


AVERAGE L= 1.188 B=. 21175 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.6639E 06 0.5964E 06 0.5050E 06 0.3904E 06 0.3108E 06 


AVERAGE L = 1.188 B=. 25851 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.8776E 06 0.7972E 06 0.6687E 06 0.4846E C6 0.3495E 06 


AVERAGE L= 1.188 B=. 25853 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.8554E 09 0.6005E 09 0.3943E 09 0.2397E 09 0.1695E 09 


AVERAGE L = 1.188 B=. 24959 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.7733E 06 0.6856E 06 0.5550E 06 0.3807E 06 0.2603E 06 


AVERAGE L = 1.188 B=. 24954 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SGUTHERN HEMISPHERE 

0.U75E 09 0.9275E 08 0.7004E 08 0.5030E 08 0.3995E 08 


AVERAGE L= 1.188 8= . 23938 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NURTHERN HEMISPHERE 

0.6690E 06 0.5769E 06 0.4486E 06 0.2890E 06 0.1862E 06 


AVERAGE L= 1.188 B=. 23934 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.1212E 08 0.1095E 08 0.9695E 07 0.83836 07 0.7571E 07 
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AVERAGE L = 1.L88 B=. 21981 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.5056E 06 0.4133E 06 0.2976E 06 0.1702E 06 0.9799E 05 


AVERAGE L = 1.188 B=. 21981 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.8335E 06 0.7838E 06 0.7133E 06 0.6223E 06 0.5579E 06 


AVERAGE L= 1.188 B=. 20482 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.4087E 06 0.3211E 06 0.2185E 06 0.1146E 06 0.6086E 05 


AVERAGE L = 1.188 B = . 20481 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.5597E 06 0.4837E 06 0.3855E 06 0.2689E 06 0.1932E 06 


AVERAGE L= 1.200 B=. 24991 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.7400E 06 0.6498E 06 0.5184E 06 0.3469E 06 0.2312E 06 


AVERAGE 1= 1.200 B=. 24992 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

C.1125E 09 0.8896E 08 0.6733E 08 0.4849E 08 0.3861E 08 


AVERAGE L = 1.200 B = . 23952 
AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.6373E 06 0.5433E 06 0.4160E 06 0.2608E 06 0. I634E 06 


AVERAGE L = 1.200 B=. 23948 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.I122E 08 0.1018E 08 0.9049E 07 0.7866E 07 0.7131E 07 


AVERAGE L = 1.200 B=. 22940 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.5506E 06 0.4569E 06 0.3357E 06 0.1976E 06 0.1166E C6 


AVERAGE L = 1.200 B=. 22938 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

C • 1 666 E 0 7 0.I638E 07 0.I594E 07 0.1536E 07 0.1496E 07 


AVERAGE L = 1.200 B=. 21979 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.4790E 06 0.3869E 06 0.2735E 06 0.I516E 06 0.8459E 05 
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AVERAGE L = 1.200 B=. 21980 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.8144E 06 0.7642E 06 0.6936E 06 0.6031c 06 0.5390E 06 


AVERAGE L= 1.200 B=. 20970 

AVERAGE NUMBER DENSITIES FOR DIFF S 
NORTHERN HEMISPHERE 

0.4136E 06 0.3249E 06 0.2207E 06 0.U51E 06 0.6061E 05 


AVERAGE L= 1.200 B=. 20971 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SOUTHERN HEMISPHERE 

0.6189E 06 0.5489E 06 0.4559E 06 0.3416E 06 0.2637E 06 


AVERAGE L= 1.200 B=. 18974 

AVERAGE NUMBER DENSITIES FOR OIFF S 
NORTHERN HEMISPHERE 

0.3116E 06 0.2325E 06 0.1468E 06 0.6659E 05 0.3260E 05 


AVERAGE L- 1.200 3=. 18975 

AVERAGE NUMBER DENSITIES FOR DIFF S 
SGUTHERN HEMISPHERE 

0.381BE 06 0.3026E 06 0.2112E 06 0.U79E 06 0.6790E 05 
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2 . Tape 10 


a. Output card description 


Columns Mode Quantity Units Description 


B, L Card 

1-8 

F 

EL 

earth radii averaged magnetic field line 


9-16 

F 

B 

gauss magnetic induction 


72-80 



constituent name 

Density 

1-12 

E 

EN1 

atoms/cm 3 density for flux model 1 

Card 

13-24 

E 

EN2 

atoms/cm 3 den " 2 


25-36 

E 

EN3 

tl H H g 


37-48 

E 

EN4 

ft ft IT ^ 


49-60 

E 

EN5 

ft ft tl 5 


72-80 

— 


constituent name 
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GENERAL PURPOSE DATA SHEET 
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GSFC FORM 541-1 (July - 40) 










H. Running Time 



An allowance of 8 seconds for each B-L line will give a close estimate of 
the running time for this program. 
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IV. LAMBDA PUNCH 


A. Introduction 

This program takes the cards from the longitudinal averaging processor 
and punches the latitude in degrees on every B-L card. This value is necessary 
in the execution of the "bounce" average calculation. 

Computing the latitude is done by the method of false position using the 
equation. 


M V4 - 3 cos 2 \ 
r e L 3 cos 6 K 


where M is the earth’s magnetic dipole moment, r & is the radius of the earth 
and k is the latitude. B is traded to the right hand side of equation (1) and the 
resulting function is evaluated for zero at a fixed B and L: 


BFUNF (B, L, k) = — tlf 3 COS — - B 
r 3 L 3 cos 6 X. 

e 

B. Mnemonics 

Quantity Description Units 

IT counter to eliminate repetition in 

computation 

LT L counter 

KT counter to determine whether a partic- 

ular L is the first L considered 

RTD conversion factor - radians to degrees 

DTR " " - degrees to radians 
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Quantity 


Description 


Units 


RE 

EM 

NPT 

BL(J) 

B(I) 

SSI (I),..., SS5(I) 
RATIO(I) 

BCOM(I) 

XLW1 

XLW2 

FXLW1 

FXLW2 

XLAMP(I) 

TERM 1 

TBLAM(J, I) 

TEST(K) 

CKT(K) 

SCKT(J) 


radius of the earth (6378.16 x 10 s ) 

the earth's magnetic dipole moment 
(8.1 x 10 25 gauss cm 3 ) divided by the 
radius of the earth cubed 

count of the number of B's to a given L 

J th magnetic field line 

I th magnetic induction for a given L 

longitudinally averaged number densi- 
ties for the five flux models and the 
I th B. 

ratio of the upper limit on the latitude 
XLW2 to the lower limit on the latitude 
XLW1 for the I th B. 

B computed for the final latitude 
XLAMP(I) at B(I) 

lower limit on the latitude 

upper limit on the latitude 

BFUNF computed for XLW1 

" » ” XLW2 

final selected latitude for B(I) 


BFUNF computed for XLW1 

latitude for the I th B of the J th L 

temporary storage of XLAMP(J) 

the absolute value of TEST(K)/TEST(K- 

temporary storage of final CKT(K) 


cm 

gauss 

earth radii 
gauss 

atoms /cm 3 


gauss 

radians 

H 

radians/ 

degrees 

degrees 

radians 
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M /4-3cos*\ 
~r^ cos® K 


C. Flow Chart 


LT = LT + 


Input and output counter NPT 
(Tapes 2, 5) 


Input BL(LT) with all the NPT B's and 
appropriate densities all with 
constituent name 



Setup constant 
EL = BL(LT) 


BC0M(J) =0 
for J = 1, NPT 


tj = f(B(J) t EL.A.j) 
ti (A.« - \i ) 

' (J) = V tl -f(B(J), EL,\ 2 )‘ 
TEST{K) = \(J) 



\ 
, 501 


mj) 


sign as \(J)? 



, 500 


\(J) 



^f(B(J),EL,\p:=0^>M = ^ 


f(B(J), EL, k 2 ) = 0? '> JL * MJ) = 


M\ m - 3 cos 2 \(J) 


Crlj L^ cos® ; 


K = K + 1 



/° 307 

CK + (K) = 

TEST (K) 1 
TEST (K - 1) j 


C 




I no 

mA 40 / 


j = j + i 


CKT(K) - 1.0001 


Cxi 

II 

*< 

yes /Does f(B(J), EL, \ 7 ) haviN™ s 


X 

the same sign as / 


J 


' Output BL{LT) with all the > 
NPT B’s latitudes and densities 
with constituent names 


Store latitudes for 
BL(LT) 


Output counter NPT 
(Tape 3) 


Output BL(LT) with all the ' 
NPT B's, latitudes, BCOM’s 
and SCKT’s 

V (Tape 3) J 



CD 




























LAMPLM PLNCH 

DJMLNSlOh ttl i'i) » TiiLAM l ?5 »^S») tSSMt b) »SS2 U •> > , SS3 US ) »;>$A (, i> ) 
lS^(2 r j) »cCGM25) 1 1 L b K 20 » »CK T i 2C ) , jLK T ( ) , XL Af'P I ? > ) 
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TfcST (a H ( J) 

IK XL /if p( J J-S I o "J F { X L A y P ( J ) , A L ^ 1 ) ) 5U 
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caTckml formula nipheks with cchrl spcnc ing internal flurula numlfhs a\u octal locations 
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E. Restrictions 


Before the first card of a new L line section a counter card must be inserted 
in order to inform the computer of the number of B's connected with that L. 

F. Input 

There are three different cards in the input deck for this program. A coun- 
ter card is followed by the designated number of card pairs for each L line of 
every constituent. The card pairs consist of a B-L card and a density card in 
that order. These are the cards output from the longitudinal averaging processor. 

1. (Tape 2) 

CARD DESCRIPTION 



Columns 

Mode 

Quantity 

Units Description 

Counter 

1-2 

I 

NPT 

B-L, density card pair counter 

Card 

B-L 

1-8 

F 

EL 

earth radii magnetic field line 

Card 

9-16 

F 

B 

gauss magnetic induction 


73-80 

- 

- 

constituent name 

Density 

1-12 

E 

SSI 

atoms/cm 3 density for flux model 1 

Card 


13-24 

E 

SS2 

ft ! 1 ft ft ft £ 


25-36 

E 

SS3 

ft ft ! T ft tt g 


37-48 

E 

SS4 

ft ft ft ft tt ^ 


49-60 

E 

SS5 

ft tt tt tt tt Cj 


73-80 

— 

— 

constituent name 
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GSFC FORM 541-1 (July • 60) 






G. Output 


Output for this program occurs on logical tapes 3 and 5. Tape 3 contains 
the results of the latitude calculation. The counter NPT is printed along with 
L, B, K, BCOMP and the ratio RATIO. BCOMP is the new B computed by in- 
serting the values of L and \ into equation (1). It should compare with the initial 
B. RATIO is the ratio of A. to a former value of the latitude computed in the 
intermediate steps. It is the same as SCKT (refer to the mnemonics listing). 
RATIO should be within .0001 of unity. 

Tape 5 contains the same data that was input on tape 2 with the exception 
that latitudes were added. This tape is punched and used as input for the "bounce” 
average calculation. 
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1. Tape 3 Sample 


NPT 

4 


L 

B 

LAMBDA 

B COMP. 

RATIO 


1. 14200 

0.23690 

9.636 

0. 2389C0CE-00 

0. 1000009E 

01 

1. 14200 

0. 22738 

7.740 

0.2273800E— 00 

0. 1000034E 

01 

1. 14200 

0.23465 

9.129 

0.2346499E— 00 

0. 1C00039E 

01 

1.14200 

0.21968 

3.667 

0. 2196800E-00 

0. 1000023E 

01 
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L 

fi 

LAMBDA 

B CUMP. 

RATIO 


1.17005 

0.23991 

12.444 

0.2399 09 9E-00 

0. 1000030E 

01 

1.17005 

0.23467 

11.752 

0.2346700E— 00 

0. 1000012E 

01 

1.1 7005 

0.22956 

11.022 

0. 229560 Oh— 00 

0. 1000005E 

01 

1.17005 

0.21976 

9.420 

0. 2197599 E— 00 

0. 1CC0056E 

01 

1.17005 

0.20973 

7.346 

0. 2097300E-00 

0. 1G00005E 

01 

1.17005 

0.20474 

6.014 

0.2Q47400h-00 

0. 1000015E 

01 

1.17005 

0.201 37 

4.894 

0 .201 3 700L-00 

0. 1000025E 

01 
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2. Tape 5 


a. Output card description 

Columns Mode Quantity Units Description 


Counter 

1-2 

I 

NPT 

- 

B-L-A. and density card pair 

Card 





counter 

B-L-\ 

1-8 

F 

EL 

earth radii 

magnetic field line 

Card 

9-16 

F 

B 

gauss 

magnetic induction 


17-24 

F 

XLAM 

degrees 

latitude 


73-80 

- 

- 

- 

constituent name 

Density 

1-12 

E 

SSI 

atoms/ cm 3 

density for flux model 1 

Card 

13-24 

E 

SS2 

H 

H 11 H M 2 


25-36 

E 

SS3 

! T 

11 IT 11 11 g 


37-48 

E 

SS4 

M 

11 It IT H ^ 


49-60 

E 

SS5 

11 

it ii H ii 5 


73-80 

- 

_ 

- 

constituent name 
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Sam P le GENERAL PURPOSE DATA SHEET 
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GSFC FORM 54 M (July - 60) 



H. Running Time 


This program will take close to a minute for every two L lines evaluated. 
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V. "BOUNCE" 


AVERAGE CALCULATION 


A. Introduction 

The "bounce" average of the number density is defined (S 5 by the equation 



s 


In other words, for a given mirror point k and a given field line L, the "bounce" 
average of the number density is the average number of atoms/cm 3 that a parti- 
cle encounters while spiraling about a field line from the northern to the southern 
mirror points (see Figure 5). The earth is assumed to be a dipole. This gives 
symmetry to the magnetic field, permitting the integrals in equation (1) to be 
evaluated over a fourth of a complete oscillation (see page 189, reference 1, 
for a discussion of the motion of trapped particles in a magnetic field.) 

The output from the lambda punch becomes input to this program. Results 
are printed and punched for use in later programs. Subroutine TABLE is in- 
cluded to interpolate (extrapolate) the density table to supply the appropriate 
density at any latitude specified by the main program. 

B. Equations 

To calculate equation (1) we project the element of arc length ds onto the 
field line L. Next we eliminate L from the equation and use the dipole 
representation 


B --^-(1 + 3 sin 2 \) 


1/2 


( 2 ) 


to express B in terms of k . This gives the "bounce" average p weighted over 
latitude for a given field line as 
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RHOAV 


( 3 ) 


r° p(k) A(\) d\ 

! 0 

f"° A(\) d\ 

Jo 

where the "weighing" factor is 

A , x „ cos 4 \ 1^4 - 3 cos 2 \ 

A(\) = - — - = 

V cos 6 \ y 4 _ 3 cos 2 \ 0 - cos 6 \ 0 V4 - 3 cos 2 A 

This factor allows for the fact that particles spiraling about the field line stay 
longer at the mirror latitudes k Q (see page 189, reference 1, for a complete 
treatment of the derivation). Figure 6 is a plot of A(\) versus k for different 
mirror latitudes. Since A(k) becomes undefined at k Q we divided equation (3) 
into two cases: 


RHOAV = 


f 


V 2h 


p(\)A(\)d\ + 




k -2h 
IL 


p(k)Pi(k)dk 


J *k -2h r k 

0 A(A.)d\ +1 A(\)d\ 


Next we designate the square of the denominator of A(A) by q(A). Expanding 
about X. Q by Taylor's series gives us 


g (\) = (k- k 0 ) 2^ — tj- 2 - (^-V 1 

i = 1 

where q(i ) (\ Q ) is the i th derivative of q evaluated at To approximate the 
infinite series we define 


3 

S (k) =SLAMF (k, k 0 ) = - 

i = l 
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where the minus sign is inserted in order to avoid the square root of a negative 
number. In the program we let S1LAMF(\, k Q ), S2LAMF(\, k Q ) and S3LAMF 
(k, k 0 ) represent the first, second and third terms of SLAMF(\, k 0 ) respectively 
and we define 


WLAMF (k, k Q ) = 


cos 4 A. Y± - 3 cos 2 \ 
VSLAMF(\, \ Q ) 


This gives 



p(\)A(\)d\ 


r 

k “2 h 


/o(\)WLAMF(\. \ 0 ) 


and 


r\> r\ 

I A(\)d K = 

\~2h k 0 - 2 h 


WLAMF (k, \ Q ) 


which are linearly approximated by the equations 


(4) 


(5) 


= 


16 


ADDNUM = y2h|j p (\ o )WLAMF(\ 0 , \) + ~ p(\ 0 -h) WLAMF (k Q - h, k Q ) 


+ P(k Q - 2h) WLAMF (A- 0 - 2h, \ Q )j 


and 


ADDDEN = V2h — 
.5 


WLAMF (\ Q , \ Q ) + — WLAMF (A 0 _ h, k Q ) + ~ WLAMF(\ Q - 2h, k Q ) 
15 
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respectively. The final step is to integrate by Simpson's rule (6) and add on 
ADDNUM and ADDDEN to get the correct density p . 

Additional equations used in this program are listed below: 

AAAF(X) = V 4 - 3 cos 2 \ 

BBBF(A) = cos 6 \ 

MONEF(X) = 3cos 4 \s i n 2\ the first derivative of BBBF(^) used in computing 

S1LAMF(A, K 0 ) 

MTWOF(X) = (cos k sin 2k) 2 - 6 cos 6 \ the second derivative of BBBF(A) 

used in computing S2LAMF(A-, ) 

MTHRF(\) = 48 cos 4 \ sin 2\ - 15 sin 3 2\ the third derivative of BBBF(^), 

used in computing S3LAMF(A.,A 0 ) 

ETAF(iV) = 10 cos 2X. - 3 sin 2 2k - 6 cos 2 2k a factor which appears in 

the second and third deriv- 
atives of AAAF(\) 

NONEF(A.) = 3 s in — the first derivative of AAAF(X) used in computing 
2AAAF(\) S1LAMF( ^ > Ao) _ 

NTWOF(X) = 3 ETAF(\) the second derivative of AAAF(\) used in 

- o cos A.) computing S2LAMF(^, k Q ). 


NTHRF(A.) 


9 ETAF(X.) sin 2k 

3 -20 sin k - 6 sin 4k - 4-3 cos 2 \ 

^ (4-3 cos 2 \) 3/2 


the third derivative 
of AAAF(\), used in 
computing S3LAMF 
(k, k 0 ) 


WGTF(\, k Q ) the designation for A(\) when integrating by Simpson's rule 
ALAM the value of the denominator in equation (3) but without equation (5) 
RHOAS the value of the numerator in equation (3) but without equation (4) 
BSUBO B Q computed by equation (2) for k Q . 
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C. Mnemonics 


Quantity 

Description 

Units 

DTR 

Conversion factor - degrees to radians 

- 

RTD 

Conversion factor - radians to degrees 

- 

RE 

radius of the earth, 6378.165 x 10 s cm. 

cm. 

EM 

earth's magnetic dipole moment, 8.1 x 10 25 gauss 
cm 3 , divided by the radius of the earth cubed. 

gauss 

HVAR 

number of partitions in Simpson's rule 

- 

NPT 

counter of the B's for a given L 

- 

XLAB1 

first six letters of constituent identification name 

- 

XLAB2 

remaining two letters of constituent identification 
name 

- 

EL 

magnetic field line L 

earth radii 

B(I) 

I th magnetic induction for a given L 

gauss 

SS(I,J) 

(a) the first time this designation is used it 
represents the longitudinally averaged density 
for the I th solar flux model and the J th magnetic 
induction 

atoms/cm 3 


(b) the second time this designation is used it 
represents the log of the densities in (a) 

- 

LCNT 

line counter for output on tape 3 

- 

XMIN 

storage of the smallest available latitude 

degrees 

XLAMW(K) 

latitude in radians, sorted in decreasing order 

radians 

XLAMRF 

storage of maximum value of the latitude 

radians 
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Quantity 

Description 

Units 

XLAMO 

mirror latitude Q 

radians 

RHOLOfl) 

density for I th flux model at mirror latitude 

atoms/ cm' 

HS 

integration interval h for Simpson's rule 

radians 

RHOAV 

"bounce" average density 

atoms/ cnr 

RATRHO 

ratio of density at mirror latitude to density 
RHOAV 

- 

RHOLOP 

density at mirror latitude 

atoms/cm‘ 

DENB(K) 

temporary storage of "bounce" average density 
at K th flux model for tape 5 punch 

atoms/ cm- 
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D. Flow Chart 

1. Main Program 


q(i) (\) is the i* 11 derivative of 
cos6 k i '4 - 3 cos2 ajj 


Read to 

counter NPT ant) 
constituent name 
(Tape 21 


Take mirror latitude \()to be that 
latitude in radians, which, when 
converted to degrees, is equal to 
the greatest integer function 
of largest input latitude 


Call Table for /O(\o)at five 
flu* models. This employs the 
stored logs f« log interpolation 
or extrapolation 


Compute ADDDEN 

f >k 0 



A(A)dA 



fN) 

cos4 A.r’4-3 cos? k 

= 

. dA 

\-0-2h 

V'(\ Q -\)S(\) 


find 2,1 A(\)d\ 

Jo 

by Simpson’s Rule 



Find Brj 

M 

A - 3 C05 2 Ajj 

T e 

co^ A„ 


Output L. k .{j. Bfi. 
and constituent name 
i Tapes 3 and 5i 


Read in the one L with its 
B’s and latitudes and five 
densities per B and 
constituent name 
(Tape 2) . 


initiation | IS/ 

interval 


/K 198 HVAR = HVAR + 2 
/ \ y" 199 


LCNT = LCNT r 1 

I I 


Cali Table for ekAfr-ht 
and PiAg-ZIU 
for model K. 


Pick ainimun latitude 
from those available and 
subscript it K 


j Convert this latitude 
to radians and store 


Are S(\ 0 ), S(Af)- h) or 
v S(\ 0 -2h) < 0? 


rNt 

=J p(>OA(.V 


Compute AD DRUM 

(>OA(M ^ 


'-0-2h 

f^o elk) cos* A A -l cos2 A. 


AjV2h V, S(>. 

for model K. 


Match up density i 

j subscripts with new k - K - 1 
; latitude subscript K } 


Begin a new page 
(Tape 31 


yes ^ Xls (LCNT -tSP ^ no 


find I a ; u A' \ .(j-. by 

J C 

Simpson's Rule for model K. 

I 

Compile RHONUM 
r'- 0 - 2 h 

Di>:, Ar \ > d A + ADDNUM 
*n for model K. 


RHONUM 
P RHODEN 
for model K. 


Output constituent name and p 
for the five flu* models 
for punching 
(Tape 51 - 


' Output K, c, o(Aq), - — — 

RHONUM, RHODEN, and third" term of 
S t A) for model K 

V. (Tape 3) 
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START 



RETURN 
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F. Input 


The input to this program are simply the cards output from the lambda 
punch. Each constituent can be run separately or they can all be run at the same 
time. Any number of L lines may be run also. A counter card must state the 
number of points to each L line which is run. Counter cards are included in the 
output from the lambda punch but they must be corrected if the number of points 
per L line are changed before running this program. Input occurs on tape 2. 

1. Input card description 


Columns Mode 


Counter 

Card 

1-2 

I 

B-L- 

1-8 

F 

Card 

9-16 

F 


17-24 

F 


73-80 

- 

Density 

1-12 

E 

Card 

13-24 

E 


25-36 

E 


37-48 

E 


49-60 

E 


73-80 

- 


Quantity 

Units 

NPT 




EL 

earth radii 

B(I) 

gauss 

XLAM(I) 

degrees 

SS(1,I) 

atoms/ cm 3 

SS(2,I) 

T 1 

SS(3,I) 

U 

SS(4,I) 

M 

SS(5 ,1) 

TT 


Description 

B-L-V, density card pair 
counter 

magnetic field line 
magnetic induction 
latitude 

constituent name 

density for flux model 1 
" " 2 

" " 3 

" *' 4 

" " 5 

constituent name 
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GSPC FORM 54 1-1 (July - 60) 



G. Output 


Output from this program appears on logical tapes 3 and 5. Tape 3 gives re- 
sults which are to be printed. Each page is devoted to a single L line for a 
specific constituent. Of the eleven columns which cross the page, the first con- 
tains the constituent name. This is followed by the value of L, the mirror lati- 
tude /v Q , the corresponding value of B, the flux model number, the "bounce” 
average number density, the longitudinally averaged number density, the ratio 
of the longitudinally averaged number density to the "bounce" averaged number 
density, the value of f 0 0 p(k)A(k) dA. , the value of f 0 0 A(\)d\ and the third 
term of SLAMF( A. 0 + 2h, A_ 0 ). All these additional values are included as an 
intermediate check on the process. 

Tape 5 contains the data to be punched. The counter card has been elimi- 
nated and all that remains are the B-L-iV card and the "bounce" average density 
card with constituent names punched to the right of each card. 
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2. Tape 5 


a. Output card description 



Columns 

Mode 

Quantity 

Units 

Description 

B-L- \ 

1-8 

F 

EL 

earth radii 

magnetic field line 

Card 

9-16 

F 

BSUBO 

gauss 

magnetic induction 


17-24 

F 

XLAMP 

degrees 

latitude 


73-80 

- 

- 

- 

constituent name 

Density 

1-12 

E 

DENB(l) 

atoms/ cm 3 

density for flux model 1 

Card 

13-24 

E 

DENB(2) 

M 

" " 2 


25-36 

E 

DENB(3) 

M 

" ” 3 


37-48 

E 

DENB(4) 

Tt 

n it 4 


49-60 

E 

DENB(5) 

IT 

ii ii 5 


73-80 

- 

- 

- 

constituent name 
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H. Running Time 


This program will take about twenty minutes for every two L lines if five 
constituents are run together. 
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VI. R, I CALCULATION AND FLUX ELIMINATION 


A. Introduction 

This program computes the atmospheric scale factor R and the atmospheric 
loss parameter 1 . The scale factor is used in the next program in order to re- 
late the energy loss of the atmosphere with the measured energy loss data. (7) 

The atmospheric loss parameter appears in the calculation of the proton loss 
term of the conservation equation in the next section. Solar flux S is eliminated 
from both R and 2 by subroutine ELIM. This subroutine logrithmically inter- 
polates Figure 9 to yield R and 2 as functions of time rather than solar flux S. 

B, L , and k and the "bounce" averaged number densities are input to this program. 
B, L , k , R and 2 are output for use in the conservation equation calculation. 


B. Equations 

For a given B and solar flux model, the "bounce" average number densities 
of the five constituents are put together to form an average number of equivalent 
oxygen atoms/cm 3 by the equation OXY = RHO/8 where 

RHO = 14n (N2 > + 8n(°> + 2b< He > + 16S<°2) + n « 

and n (J ) j s the "bounce" average number density for the I th constituent. 

The scale factor R is given by the equation 


D ..x _ (OXYGEN ATOMS /CM 3 ) ATMOS 

Kylw, D, t) - 

(OXYGEN ATOMS /CM 3 )NTP 

where (oxygen atoms/cm 3 ) ATMOS is OXY and (oxygen atoms/cm 3 ) NTP comes 
from the following relationship of an ideal gas: 


22414cm 3 /Kmole = .6 024 9 x 10 24 atoms/Kmole 

or (oxygen atoms/cm 3 ) NTP = 2.69 x 10 19 . Figure 7 is an example of the output 
from this program. It shows the time dependence of the atmosphere in terms of 
the scale factor R . 
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Sigma is given by the equation 




.-(He) + 



+ 2n (0 2 ) 
8 



cr(0) atoms/cm. 


where ^(He)and ^(0) are the interaction cross sections of helium and oxygen 
respectively. Figure 8 again illustrates output from this program. It shows 
log e I as a function of time for various values of B at anL of 1.25 earth radii. 


C. Mnemonics 

Quantity 

NPT 

C 

SIGHE 
SI GO 
EL 
B(N) 

ALATO(N) 

HE1(N),. . ., HE5(N) 

01 (N), . . ., 05(N) 

021 (N), .. ., 02 5 (N) 

AN21(N), . . ., AN25(N) 


Description 


counter of the number of B’s to a 
given L 

(oxygen atoms/cm 3 ) NTP = 2.69 x 10 19 

cr(He) = .143 x 10- 24 

o-(O) = .36 x 10“ 24 

magnetic field line L 

magnetic induction B for a 
given L line 

latitude corresponding to B(N) 

helium "bounce" averaged densities 
for B(N) and the five flux models 

oxygen "bounce" averaged densities 
for B(N) and the five flux models 

molecular oxygen "bounce" averaged 
densities for B(N) and the five flux 
models 

nitrogen "bounce" averaged densities 
for B(N) and the five flux models 


Units 


atoms /cm 3 
cm 2 

ft 

earth radii 

gauss 

degrees 

atoms /cm 3 

atoms /cm 3 

atoms/cm 3 

atoms/cm 3 
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Quantity 

Description 

Units 

H1(N), . . H5(N) 

hydrogen "bounce" averaged densities 
for B(N) and the five flux models 

atoms /cm 3 

RH01(N), . . .,RH05(N) 

RHOfor B(N) and the five flux models 
(see equation on page 105) 

atoms/cm 3 

OXYl(N), .... OXY5(N) 

OXY for B(N)and the five flux models 
(see equation on page 105) 

atoms/cm 3 

SIG1(N), . . SIG5(N) 

atmospheric loss parameter 2 for 
B(N) and the five flux models 

atoms /cm 

RATI (N), . RAT5(N) 

scale factor E for B(N) and the five 
flux models 

- 

LCNT 

counter to notify subroutine whether 
it is working with R or 2 

- 

ARG1 

temporary storage of RATI or SIG1 
depending on LCNT 

depends on 
LCNT 

ARG2 

temporary storage of RAT2 or SIG2 
depending on LCNT 

depends on 
LCNT 

ARG3 

temporary storage ofRAT3or SIG3 
depending on LCNT 

depends on 
LCNT 

ARG4 

temporary storage ofRAT4or SIG4 
depending on LCNT 

depends on 
LCNT 

ARG5 

temporary storage ofRAT5or SIG5 
depending on LCNT 

depends on 
LCNT 
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2. Subroutine Elim. 


START 
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RH04(*. ) = L4.*AN?MM+8.»CMi\)+2**HL<»(M + l6.«C?MN)+H*t(iN) 
Rh05(\) = lA.*AN25tM+8.*05tN)+2.«HLt>(N» + lfc,»C^i)<N)+HblN) 
CCNT I\tf: 
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ANS = Y C J ) 

(iO TO 1CCO 
H 1 = S ( J ) 
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STORAGE LOCATIONS FCR VARIABLES APPEARING IN CCNKCN STATEMENTS 
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F. Input 


Input to this program are the "bounce" average density cards together with 
the B-L- A. cards output from the last program. The cards are rearranged by 
L line rather than by constituent as they have been arranged in previous programs. 
A counter card must again be inserted in order to inform the computer of the 
number of card pairs per constituent for a , given L- line. All input is on tape 2. 

For a given L-line the constituents must be in the following order: 

(1) helium, 

(2) oxygen, 

(3) molecular oxygen, 

(4) nitrogen, 

(5) hydrogen. 


1. Input Card Description 



Columns 

Mode 

Quantity 

Units 

Description 

Counter 

Card 

1-2 

I 

NPT 

- 

B-L- A. and density card 
pair counter 

B-L- A. 

1-8 

F 

EL 

earth radii 

magnetic field line 

Card 

9-16 

F 

B 

gauss 

magnetic induction 


17-24 

F 

ALATO 

degrees 

latitude 


73-80 

- 

- 

- 

constituent name 

Density 

1-12 

E 

* 

atoms/cm 3 

density for flux model 1 

Card 

13-24 

E 

* 

ft 

density for flux model 2 


25-36 

E 

* 

ft 

density for flux model 3 


37-48 

E 

* 

ft 

density for flux model 4 


49-60 

E 

* 

It 

density for flux model 5 


73-80 

_ 

_ 

- 

constituent name 


All of the five constituents use the same card format here. Therefore the quantity on a given 
density card (HE, 0, 02 N2 or H) will depend on the constituent name listed in columns 73-80 
of that same card. 
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G. Output 


Output for this program occurs on tapes 3 and 5. Tape 3 contains two 
groups of data for each L. The first group lists the values of R , I andOXY as 
functions of position and flux model. The second group contains R, log e R, 2, 
log e 2 and solar flux as functions of position and time. 

On tape 5 is the data to be punched for use in the conservation equation cal- 
culation. Each B-L-X card is followed by two R cards which are followed in turn 
by two 2 cards. 
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1. Tape 3 Sample 



8 = 750 

s=2'<u 

S = 1 5 0 

S=100 

$=70 



L= 1.14203 

o 

v> 

T' 

■\J 

• 

II 

Xj 

LAT= 6. 

0-73 



'“'XY 

9.34/92 11 

n. 39 7 f > r. i ] o 

.0624 2 11 

0. 8475E 1 1 

0 • 1 0 A 6 i : 

12 

A T I U 

0 .2 0 3 ( b - 0 o 

0.2221E-00 0 

• 2637E-0B 

y. 3131E-0G 

0. 3 8 892 

-OB 

SIGMA 

0.30 ilu-14 


. 3 7 602-14 

0.467 22-14 

0. 37:371: 

“14 



L = 1.14 20 

U= .1.22} 94 

L A T - 6. 

9 7 3 


GXY 

0.3230t 07 

D . 2 3 4 6 C U 9 

0 • 1 b 3 3 2 o y 

0.1138E 00 

0. 74320 0* 

RATIO 

1 2 0 8 2-19 

0.9462E-1 l 

0. 6822c- L l 

0.42312-1 1 

3.27682-11 

SIGMA 

0. 1 763E-10 

0. 1 J7 Jfc-16 

0.9906E-17 

0 . 6 1 0 1 E - 1 7 

0. 39348-1 7 



L - 1.14200 

0= 0.21663 

L A T = 4. 

>73 


□ XY 

0 . 1 6 3 1 l 0 H 

0.9)6';,- 07 

0.32782 07 

O. 21 41 h 07 

0.1007c 07 

RATIO 

0 . 6 0 6 5 1 - 1 2 

0. 1 7 0 0 E - 1 2 

0. 19622-12 

0 ./ 9 5 8 2 - 1 3 

0.3744c-l 3 

SIGMA 

0.8 4 Ooc -18 

0. *j144<>1E 

0 . ? 7 4 7 L - i 6 

0. 1 1 4 4 1 - 1 8 

0 . 3 6 3 3 2 - 1 9 



L = 1.14200 

U- 0.2121 

6 L A r = 2 . < 

4 7 3 


0XY 

0.1 7 TIE 0 7 

0.6907c 06 

0. 3902c Oo 

0 .12 9 7 E 06 

0.34242 03 

RATIO 

0 . b 4 3 3 E - 1 3 

0.334-5 c -13 

0. 1 4 b U h - 1 3 

0 . 4 81 6 Q - 1 4 

0 * 2 0 1 m r - 1 4 

SIGMA 

0 • 10 2 3 1 - 1 c 

j . ■-> ) 0 y E “ 1 0 

0 . 3 1 7 3 c - 1 9 

3.14252-19 

0. 7 H 7 0c“2^ 



L= 1.14230 

L - 0.2096 

5 L AT = 0. 

y 73 


□ XY 

0.3371E 06 

0.1677L Oo 

0 • 7 3 9 8 t 0 3 

0.2 98 IE 0 3 

0.1367c 05 

RATIO 

0. 126 IE-13 

0 . 6 2 3 2 E - 1 4 

0 . 2 0 2 4 E - 1 4 

0. 11U8E-14 

0. 3 7 8 v c - 1 3 

SIGMA 

0. 3246E-1 ) 

(J.2i^32-l7 

0.1314h-l y 

0.6010 2-20 

0. 392 72-20 
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LAT = 

a. 9730 L= 

1.1420 B= 

0.2338 


"I ME 

SOLAR 

LUG UF 

RATIO 

YRS) 

FLUX 

RAT I U 



0 

70.00 

-0.1 9365F 

02 

0. 38892E-08 

1 

73.00 

-0. 19400E 

02 

0. 37550E-08 

2 

130.00 

-0. 19706E 

02 

0.27666E-08 

3 

230.00 

-0. 19977E 

02 

0.21 087E-08 

4 

230.00 

-0.20012L 

0? 

0. 20368E-08 

5 

220.00 

-0.1 99t»0 E 

02 

0.21 457E-0B 

6 

185.00 

-0 • 1 98855 

02 

0.23118E-08 

7 

140.00 

-0.197495 

02 

0. 264925-08 

8 

103.00 

-0. 19597E 

02 

0.30832E-06 

9 

90.00 

-0. 1 95 J 35 

02 

0.33 798E-08 

10 

75.00 

-0. 19400E 

02 

0.37550t-06 

11 

70.00 

-0. 19365c 

02 

0.388925-08 


TIME 

SOLAR 

LUG UF 


SIGMA 

(YRS) 

FLUX 

SIGMA 



0 

70.00 

-0. 32788b 

02 

0.5 7570E-14 

1 

75.00 

-0. 32823c 

02 

Q.55601E-14 

2 

130.00 

-C.33126E 

02 

0. 41067 F. -14 

3 

230.00 

-0 .333965 

02 

0.31370E-14 

4 

250.00 

-0.334305 

02 

0.303075-14 

5 

220.00 

-0.33378L 

02 

0. 319155-14 

6 

1 8 5 • 0 0 

-U. 33304c 

02 

0.343655-14 

7 

140.00 

-0.33169E 

02 

0.39339E-14 

8 

103.00 

-0. 3301 9 E 

02 

0.45726F-14 

9 

90.00 

-0. 329285 

02 

0.30088E-14 

10 

75.00 

-0. 32B23E 

02 

0.356 01E-14 

11 

70.00 

— U • 3 2 7 h 8 c 

02 

J. 375705-14 


L AT = 

6.9730 L = 

1.1420 B = 

C. 2239 


TIME 

SOLAR 

LOG OF 

RATIO 

(YRS) 

FLUX 

PAT I u 



0 

70.00 

-0.266071: 

02 

0.278325-11 

1 

73.00 

-0.?6537l 

02 

0 • 29862F-1 1 

2 

130.00 

-0. 2i> : )02t 

0 2 

0.563365-11 

3 

2 3 j . 0 0 

-0.2 r237c 

02 

0. 10957l-10 

4 

230.00 

-0.2P139E 

02 

0. 120835-10 

5 

220.00 

-0. 2S286E 

02 

0. 10434E-10 

6 

183.00 

-0.2r.482F 

02 

0.857/45-11 

7 

140.00 

-0.2r806E 

02 

0. 62006E-11 

8 

103.00 

-0.2 0 14 IE 

0? 

0.4*3815-11 

9 

90.00 

-0.2432SC 

0 2 

O'. 36806E-11 

10 

75.00 

-G.26837E 

02 

0.298625-11 

11 

70.00 

-0.2660 7 C 

02 

0.278525-11 
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"■ I ME 

sola* 

LOG OF 


SIGVA 

YRS) 

FLUX 

SIGMa 



0 

70.00 

- 0 . 4 0 0 6 4 l 

02 

0. 398446-1 7 

1 

7b. 00 

-0. 3 9 9 9 3 F 

02 

0*4277bt-17 

2 

no. oo 

-0.3‘9347fc: 

02 

0. 81603E-1 7 

3 

230.00 

-0. 3 8 6 7 5 E 

02 

0.159916-16 

4 

250.00 

-0.385766 

02 

U. 176486-16 

5 

220.0 0 

-0.36724L 

02 

0. 152216-16 

6 

18b. 00 

-0.38922b 

02 

0. 124886-16 

7 

14 0.00 

- 0 . 3 > 2 b 0 f: 

02 

0.899116-17 

8 

105. on 

—0.3959 Ob 

02 

0.64039F-17 

9 

90. OJ 

-0.39780L 

02 

0.529316-17 

10 

7b. 00 

-0.39993b 

0 2 

0.427766-17 

11 

70.00 

-0.40064c 

02 

0. 3 9 8446-17 


L AT = 

4.9730 L~ 

1.1420 6= 

0 • 2 1 6 8 


TIME 

SOLA* 

LOG uF 


RATIO 

( YRS) 

FLUX 

RATIO 



0 

7 0 . 9 0 

-0. 3091m. 

0 2 

0.374416-13 

1 

7 9.00 

-0 • 307906 

02 

0. 424556-13 

2 

130.00 

-0.296206 

02 

0. 136776-12 

3 

230.00 

-0.28329b 

0 2 

0.49 772E-12 

4 

230. 00 

-0.28 1 3 1 F 

02 

0.606496-12 

5 

220.00 

-0.2842bt 

02 

0.45089E-12 

6 

18b. 00 

-0.2 o 8 1 5 £ 

02 

0. 30591E-12 

7 

140.00 

-0. 294406 

0 2 

0. 163826-12 

8 

105.00 

-0.300726 

9? 

0.871006-13 

9 

90.00 

-0.304136 

0 2 

0.618976-13 

10 

7b. 00 

-0. 307906 

0? 

0.42455C-13 

11 

70.00 

-0.30916b 

02 

0.374 41b- 13 


TIME 

solar 

LOG OF 


SIGMA 

(YRS) 

FLUX 

SIGMA 



0 

70.00 

-0.443236 

02 

0.56327E-19 

1 

75.00 

-0. 44205b 

02 

0.63395E-19 

2 

130.00 

-0 .430896 

02 

0.193566-18 

3 

230.00 

-C. 418176 

02 

0.69067E-18 

4 

250.00 

-0. 416206 

02 

0.84064E-18 

5 

220.00 

— 0.41915E 

02 

0.62604E-18 

6 

185.00 

-0.423006 

0? 

0.426136-18 

7 

140.00 

-0.429146 

02 

0.23059E-18 

8 

105.00 

-0. 435266 

02 

0.124956-18 

9 

90.00 

-G.43850E 

02 

0.90378E-19 

10 

75.00 

-0.44205E 

02 

0.63395E-19 

11 

ro.oo 

-0.44323E 

02 

0. 56327E-19 
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LAT = 

2.9730 L= 

1.1420 6 = 

0.2121 


T I ME 

SOLAR 

LUG L)F 

RATIO 

YRS ) 

FLUX 

RATIO 



0 

70.00 

-0.33838b 

02 

0.201 6 3E -14 

1 

75.00 

-0.33692E 

02 

0.2331 lt-14 

2 

130.00 

-0.32293E 

02 

0.94462t-14 

3 

230.00 

-0. 39636E 

02 

0.495286-1 3 

4 

250.00 

-0.30374E 

02 

0.64346t-l 3 

5 

220.00 

-0.30767E 

02 

0.434526-13 

6 

185.00 

-0.31273E 

02 

0.26189F-13 

7 

140.00 

-0. 32069b 

02 

0.1l824h-i3 

8 

105.00 

-0. 32855b 

02 

0.53864F-14 

9 

90.00 

-0. 33257b 

02 

0.36028E-14 

10 

75.00 

-0. 33692E 

02 

0.233118-14 

11 

70.00 

-0. 33B3«b 

0? 

0.201636-14 


TIME 

SOLAR 

LOG OF 

SIG*4 

C YRS) 

FLUX 

SIGMA 



0 

70.00 

-C. 46290b 

02 

0 • 7 ft 7 8 3 1 - 2 0 

I 

75.00 

-0.46191b 

02 

0. 869676-20 

2 

130.00 

-0.45217E 

02 

0.23039E-1 9 

3 

230.00 

-0.4 3942E 

02 

f). 824576-19 

4 

250.00 

-0.43727E 

02 

0. 10228t-18 

5 

22C.00 

-0.4405 or 

02 

0. 740376-19 

6 

185.00 

-0.44455b 

02 

0.493816-19 

7 

140.00 

-0 • 45Gd 7E 

02 

0.27037E-19 

8 

105.00 

-0.45617E 

02 

0 • 1 5 4 4 2 E - 1 9 

9 

90.00 

-0.45395b 

02 

u . 1 1 69 BE — 1 9 

10 

7 d . 0 0 

-0. 4o 19 1C 

02 

0. R6967L-20 

11 

7 0.00 

-0.4S290E 

02 

0.787636-20 


LAT = 

0.9730 L = 

1.1420 n= 

0.2099 


TIME 

SOLAR 

LUG OF 

RATIO 

(YRS) 

FLUX 

RATIO 



0 

70.00 

- 0 • 3 3 0 3 5 l 

02 

0.578888-15 

1 

7 5.00 

-0.34977c 

02 

0. 645046-15 

2 

130.00 

-0.3 *87 3b 

0? 

0.194268-14 

3 

230.00 

-0. 3 ? 2 9 0 b 

02 

0.94762E-14 

4 

230.00 

-0.32011c 

0 2 

0. 12530F-I3 

5 

220.00 

-0. 324306 

02 

0.824096-14 

6 

185.00 

-0. 32946b 

02 

0.49151L-14 

7 

1 4 u * 0 0 

-0. 33688E 

02 

0. 23424 E; -14 

8 

} 05.00 

-0. 34343b 

0? 

0.121676-14 

9 

90.00 

-0.34653b 

02 

0.892418-1^ 

10 

7 ‘j . 0 0 

-0.34977E 

02 

0 • 64504b- 1 5 

11 

7 0 • 0 0 

-0. 35065E 

02 

0 * 5 7 8 8 8 £ - 1 5 
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^IMb 

SOLA 9 

LOG OF 


7RS ) 

FLUX 

SIGMA 


0 

70.00 

-0.4'.9bt>E 

0 2 

1 

73.00 

“0. 468930 

02 

2 

130.00 

“ U • 4 6 041 K 

0 2 

3 

230. Ou 

-U.43039L 

02 

4 

230.00 

-0.44872c 

02 

3 

2 2 0 . U U 

-0.49122L 

02 

6 

183.00 

-0.43436 1; 

0 2 

7 

140.00 

-0.45910b 

02 

8 

10 ‘>.00 

“0# 4 6 3 7 0 £ 

02 

9 

9 1 . 0 U 

“ 0 • 4 6 6 1 9 L 

0 2 

10 

7 3.00 

“0 . 403337 

0 2 

11 

7 0 . 00 

- 0 . 4 3 9 H 6 F 

0 2 



SIGMA 

0. 39268F-20 
0. 43U43L-20 
0.10104F-19 
0.27d 41L“19 
0 • 3 2 3 3 0 L - 1 9 
0.29333L”19 
0. 18399F “19 
0* 113241 -19 
0. 72732r:“20 
0.36bS5L-20 
0.4 30431- -20 
0 • 3 9 2 6 8 L - 2 0 
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2. Tape 5 


a. Output Card Description 



Columns 

Mode 

Quantity 

Units 

Description 

B-L-A. 

1-8 

F 

ALATO 

degrees 

latitude 


9-16 

F 

EL 

earth radii 

magnetic field line 


17-24 

F 

B 

gauss 

magnetic induction 

1st R 

1-10 

E 

XS(1) 

- 

R for time 0 

Card 

11-20 

E 

XS(2) 

- 

R for time 1 


21-30 

E 

XS(3) 

- 

R for time 2 


31-40 

E 

XS(4) 

- 

R for time 3 


41-50 

E 

XS(5) 

- 

R for time 4 


51-60 

E 

XS(6) 

- 

R for time 5 

2nd R 

1-10 

E 

XS(7) 

- 

R for time 6 

Card 

11-20 

E 

XS(8) 

- 

R for time 7 


21-30 

E 

XS(9) 

- 

R for time 8 


31-40 

E 

XS(10) 

- 

R for time 9 


41-50 

E 

XS(ll) 

- 

R for time 10 


51-60 

E 

XS(12) 

- 

R for time 11 

1st 2 

1-10 

E 

XS(1) 

atoms/cm 

2 for time 0 

Card 

11-20 

E 

XS(2) 

n 

2 for time 1 


21-30 

E 

XS(3) 

m 

2 for time 2 


31-40 

E 

XS(4) 

ii 

2 for time 3 


41-50 

E 

XS(5) 

ii 

2 for time 4 


51-60 

E 

XS(6) 

ii 

2 for time 5 
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Columns 

Mode 

Quantity 

Units 

Description 

2nd 2 

1-10 

E 

XS(7) 

atoms/cm 

2 for time 6 

Card 

11-20 

E 

XS(8) 


2 for time 7 


21-30 

E 

XS(9) 

t t 

2 for time 8 


31-40 

E 

XS(10) 

ft 

2 for time 9 


41-50 

E 

XS(ll) 

ft 

2 for time 10 


51-60 

E 

XS(12) 

ft 

2 for time 11 


This five card group is repeated for each new B and L . 
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H. Running Time 


The running time for this program will be about three minutes for every L 
line. 
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VH. CONSERVATION EQUATION CALCULATION 


A. Introduction 

This program studies the build-up of proton density by use of the conserva- 
tion equation (see equation (1)). Several lines of action are available. Either a 
transient or a time-averaged steady state solution can be calculated for any of 
15 desired energy levels without altering the program. 

In the transient steady-state solution the conservation equation is integrated 
as a function of time for a particular energy level until maximum and minimum 
values of the densities start repeating from one solar cycle to another. At this 
point the cycle number is recorded together with the time of the maximum and 
minimum densities. These values are printed with the maximum and minimum 
density and flux values as well as with the ratio of the maximum flux to the min- 
imum flux. In addition, when calculating the transient steady-state solution, one 
may have the time history of any one of the energy levels printed on tape 5 with 
as many points as desired. 

In the time-averaged steady-state solution the condition dN p /dt = 0 yields a 
density equation for density N p which is solved using the relative neutron source 
strength $ together with I and R from the last program. All these values are 
averaged over time for this solution of the conservation equation. 

Three subroutines assist the main program: subroutine RUNGE is used in 
integrating by the Runge Kutta technique, subroutine DERIV is used to evaluate 
dNp/dt for a given t and N p and subroutine TABLE interpolates any given table 
either logarithmically or linearly. 


B. Equations 

The form of the particle conservation equation used for the study of the 
proton population as a function of time is given by 


dN 

p 

dt 


= C n D - C N 

0 1 P 



N — 
2 p dE 



C 2 N P S 


( 1 ) 
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where — , — ( — \ $ and I are the quantities to be supplied and where: 
dx dE \dx j 

N = density 

p 

C 0 =A 0 /L 2 E B °cos 4 \ 0 
C i =A i/ £Bl 

c 2 =a / 2 

E = Energy 

A Q , . . . , A 2 = high or low energy conservation equation coefficients depending on 
whether E > 80 Mev. or E 1 80 Mev. respectively 

r> o — n tt n it m it if it fi if ft ii ti 

B 0- • • • - B 2~ 

A. = mirror latitude 

The condition dNp/dt = 0 yields the equation for the time averaged steady 
state proton density: 


A 0 $ 


l 2 e b ° 


(-) 

+ ak b > 1 

l^i) + A E B 2 Z 

cos 4 X. 


_E Bl 

\dx) 

2 dE 

\dXj 2 

0 


where $ , 2 and the scale factor R used in calculating dE/dx and d/ dE (dE/ dx ) 
are all averaged over time. 

The density N p gives a flux by the equation: 

Flux = 2C N = N v 

2 p p 

where v is the neutron velocity factor. 


132 


C. Mnemonics 


Quantity 

Description 

Units 

TIME(J) 

Abscissa of time in increments of years for 
12 years 

years 

ELOSS(J) 

dE/dk corresponding to E(J) (see Figure 11) 

Mev/ cm. 

E(J) 

energy corresponding to ELOSS(J) " 

Mev. 

CONVM 

conversion factor - months to seconds 

- 

ALO 

A 0 for E 1 80 Mev. (see page 130) 

# protons 
cm. sec. Mev. 

AL1 

A, " ” " 

cm. /sec. 

AL2 

A M IT IT 

A 2 

cm./Mev 2 /sec. 

AHO 

A 0 for E > 80 Mev. (see page 130) 

# protons 
cm. sec. Mev. 

AH1 

^ it it it 

cm. /sec. 

AH2 

A 'I ii ii 

A 2 

cm. /Mev 2 /sec. 

BLO 

B q for E £ 80 Mev. (see page 130) 

- 

BL1 

r> t! n n 

D 1 

- 

BL2 

TD 1! 11 11 

n 2 

- 

BHO 

B 0 for E > 80 Mev. (see page 130) 

- 

BH1 

g 11 11 11 

- 

BH2 

g 11 11 M 

- 

DELOSS(J) 

d(dE/dx)/dE corresponding to DE(J) 
(see Figure 12) 

cmr 1 

DE(J) 

energy corresponding to DELOSS(J) 
(see Figure 12) 

Mev. 
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Quantity 

Description 

Units 

PREL(J) 

'f' for TIME(J) (see page 130 and Figure 11) 

- 

AVPRL 

simple average of PREL(l), .... PREL(12) 

- 

ALAT 

mirror latitude \ Q 

degrees 

EL 

magnetic field lineL 

earth radii 

B 

magnetic induction B 

gauss 

CONVR 

conversion factor — degrees to radians 

- 

ALATO 

ALAT in radians 

radians 

CSLAT4 

eos 4 (ALATO) 

- 

ELEL 

L 2 

(earth radii) 2 

R(I) 

atmospheric scale factor R for TIME (I) 

- 

SIG(I) 

atmospheric loss parameter 2 for TIME (I) 

atoms/cm. 

ENO 

initial density N po for integration 

# protons/cm.3 

TSUBO 

initial time t o for integration 

months 

DTI 

integration interval for ICSUBO 

months 

DT2 

integration interval for all other cycles 

months 

TEND 

end limit on integration 

months 

IEST 

initial energy level subscript 

- 

IEDEL 

increment for energy level subscript 

- 

IEEND 

final energy level subscript 

- 

IE PR 

subscript of energy level whose time history 
is desired 
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Quantity 

NPRINT 

ICSUBO 

SUMR 

AVR 

SUMSIG 

AVSIG 

T 

TMAX 

TCK 

DT 

ICYCLE 

EMAX(J) 

EMIN(J) 

ETMAX(J) 

ETMIN(J) 

M 

ENER 

ITEST3 


Description 


a control factor for the time history print on 
tape 5. Density and flux are given for every 
NPRINT increments in time. (See Restriction 
( 1 )). 

initial cycle number corresponding to ENO and 
TSUBO 

sum of R(l), R(12) 

simple average of R(l), R(12) 

sum of SIG(l), SIG(12) 

simple average of SIG(l), SIG(12) 

time 

time of an eleven year cycle — 132 months 
time check to see if TEND has been reached 
increment in time 
cycle number 

maximum density for cycle J 
minimum density for cycle J 
time at EMAX(J) 
time at EMIN(J) 

subscript to indicate the energy level which is 
under consideration 

energy level M under consideration 

test value to see if the energy level under 
consideration is to have its time history 
printed on tape 5 


Units 


atoms/ cm 

atoms/cm 

months 

months 

months 

months 

atoms/cm 3 
atoms/ cm 3 
months 
months 

Mev. 
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Quantity 
A LOSS 
X(l), 

X(15) 

Y(l), .... 
Y(15) 

ALOSSA(J) 

DLOSSA(J) 

AO 

A1 

A2 

BO 

B1 

B2 

CO 

Cl 

C2 

DEDX 

DDEDX 


Description 

temporary storage of ELOSS(M) 


Units 


Mev./cm. 


temporary storage of DE(1), DE(15) Mev 


temporary storage of DELOSS(l), ..., 
DELOSS(15) 

(dE/dx) x R(J) for energy ENER 
(d(dE/dx)/dE) x R(J) for energy ENER 
temporary storage of ALO or AHO 

temporary storage of AL1 or AH1 

temporary storage of AL2 or AH2 

temporary storage of BLO or BHO 

temporary storage of BL1 or BH1 

temporary storage of BL2 or BH2 

conservation equation coefficient 
(see page 129) 


Mev./cm. 

cm -1 

# protons 
cm. sec. Mev. 

cm/ sec 

cm/Mev 2 sec 


# protons /cm 3 
sec Mev. 


conservation equation coefficient 
(see page 129) 

conservation equation coefficient 
(see page 129) 

temporary storage of ALOSSA(J) 
temporary storage of DLOSSA(J) 


cm/sec. Mev. 

cm/sec. Mev. 

Mev./cm. 

cm -1 
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Quantity 

Description 

Units 

DENS 

proton density N p 

# protons/cm 3 
Mev. 

FLUX 

proton flux 

# protons/cm 2 
sec. Mev. 

NP 

counter of the increments in time; used as a 
check for NPRINT 

- 

N5 

counter of the number of prints on tape 5 

- 

EN 

proton density N p 

atoms/ cm 3 Mev 

ICYM1 

the value ICYCLE-1 used in testing for 
steady-state 

- 

TEST1 

test to see if steady state has been reached 

- 

TEST2 

test to see if steady state has been reached 

- 

EPR1 

minimum density in cycle ICYM1 

# protons/cm 3 
Mev. 

FLPR1 

flux for EPR1 

# protons/cm 2 
sec. Mev. 

TPR1 

time of EPR1 

months 

EPR2 

maximum density in cycle ICYM1 

# protons/cm 3 
Mev. 

FLPR2 

flux for EPR2 

# protons/cm 2 
sec. Mev. 

TPR2 

time of EPR2 

months 

RATIO 

ratio of FLPR2 to FLPR1 

_ 
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D. Flow Chart 



L38 

















2. Subroutine RUNGE 
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3. Subroutine DERIV 
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START 
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AL0=CCi*vM 2. 56 4 E- 13 
ALl =CL'mVK» 3 .4t3fc+fc 
AL2 = CC n V ^ «7.z l 3 , 5t+t 
BL 0 = 2. 3tS 
BL 1 = . ml 3 
BL2=.4 i 1 
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LUC AT I CN S UF NAMES IN TRANSFER VECTCR 
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CALL HL«iV( TC ♦YO,CCR) 
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STCRAGc LOCATIONS FCH VARlAtiL-S APPEALING IN CCPPCN STATEMENTS 
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SUBRCLHNt TU INTERPCLATC A GIVEN TABLE 
SUBRCli I INC TABLh { XARGtVANb,LL) 

Clf'ENS IUN ) ,Y ( IS) 

CCMNflN CC.Cl,C2»X f Yf tl Ht , MR t L . ALCSSA t LLC S i» A , ^ I G 
IHLLUltllf 12 
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F. Restrictions 


(1) 132/(NPRINT x DTI) must be less than 5000 in order that T5, E5 and 
FL5 subscripts do not exceed their dimension specifications. 

(2) IEEND cannot exceed 15 since the table of energy levels contains only 
15 values (see energy level table below) 


G. Input 

Input to this program consists of five data cards and four control cards for 
every B-L line considered. The data cards are those output by the preceeding 
program. The control cards are those that handle the selection of the many 
available calculations in this program. The first control card contains the ini- 
tial density (END) and the initial time (TSUBO) for integration. The second 
card contains the integration interval for the first cycle (DTI), the integration 
interval for all other cycles (DT2) and the end limit of integration (TEND). 
Setting DTI equal to zero will give a time averaged steady state solution of the 
equation. The next control card defines the energy level(s) to be used. It gives 
the initial energy level subscript (IEST), the increment in energy level sub- 
scripts (IEDEL) and the final energy level subscript (IEEND). Energy levels 
are as follows: 


E(l) = 10 Mev 
E(2) = 25 Mev 
E(3) = 50 Mev 
E(4) = 75 Mev 
E(5) = 100 Mev 


E(6) = 125 Mev 
E(8) =150 Mev 
E(8) = 200 Mev 
E(9) = 250 Mev 
E(1 0) = 300 Mev 


E(ll) = 350 Mev. 
E(12) = 400 Mev. 
E(13) = 500 Mev. 
E(14) = 600 Mev. 
E(15) = 700 Mev. 


Note that IEST, IEEND and IEDEL must be greater than zero and that IEND 
cannot be greater than 15. The fourth control card contains three numbers (1) 
IEPR is the subscript of that energy level whose time history is to be printed 
on tape 5. If IEPR = 0 then no time history will be printed. (2) NPRINT tells 
the computer how often to print density, flux and time values on tape 5 if 
IEPR ^ 0. Values will be given for every NPRINT increments in time. How- 
ever, if 132/(NPRINTxDTl) is greater than 5000 the computer will print an 
error message and go on to the next case. This is done in order to prevent 
time, energy and flux subscripts from exceeding their dimension statement 
capacities. (3) ICSUBO is the cycle number at initial density (ENO) and initial 
time (TSUBO). ICSUBO is always greater than or equal to 1. All input occui's 
on tape 2. 
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1 . Input Card Description 



Columns 

Mode 

Quantity 

Units 

Description 

Card 1 

1-8 

F 

ALAT 

degrees 

latitude 


9-16 

F 

EL 

earth 

radii 

magnetic field line 


17-24 

F 

B 

gauss 

magnetic induction 

Card 2 

1-10 

E 

R(l) 

- 

scale factor R for TIME(l) 


11-20 

E 

R(2) 

- 

scale factor R for TIME (2) 


21-30 

E 

R(3) 

- 

scale factor R for TIME(3) 


31-40 

E 

R(4) 

- 

scale factor R for TIME(4) 


41-50 

E 

R<5) 

- 

scale factor R for TIME (5) 


51-60 

E 

R(6) 

- 

scale factor R for TIME (6) 

Card 3 

1-10 

E 

R(7) 

- 

scale factor R for TIME(7) 


11-20 

E 

R(8) 

- 

scale factor R for TIME (8) 


21-30 

E 

R(9) 

- 

scale factor R for TIME(9) 


31-40 

E 

R(1 0) 

- 

scale factor R for TLME(IO) 


41-50 

E 

R(ll) 

- 

scale factor R for TIME (11) 


51-60 

E 

R(12) 

- 

scale factor R for TIME(12) 

Card 4 

1-10 

E 

SIG(l) 

atoms/ cm 

atmospheric loss parameter 2 
for TIME(l) 


11-20 

E 

SIG(2) 

atoms/cm 

atmospheric loss parameter 2 
for TIME(2) 


21-30 

E 

SIG(3) 

atoms/cm 

atmospheric loss parameter 2 
for TIME (3) 


31-40 

E 

SIG(4) 

atoms/cm 

atmospheric loss parameter 2 
for TIME(4) 


41-50 

E 

SIG(5) 

atoms/ cm 

atmospheric loss parameter S 
for TIME(5) 


51-60 

E 

SIG(6) 

atoms/cm 

atmospheric loss parameter 2 
for TIME (6) 
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Columns 

Mode 

Quantity 

Units 

Description 

Card 5 

1-10 

E 

SIG(7) 

atoms/cm 

atmospheric loss parameter 1 
for TIME (7) 


11-20 

E 

SIG(8) 

atoms /cm 

atmospheric loss parameter £ 
for TIME(8) 


21-30 

E 

SIG(9) 

atoms/cm 

atmospheric loss parameter 2 
for TTME(9) 


31-40 

E 

SIG(10) 

atoms/cm 

atmospheric loss parameter 2 
for TIME(10) 


41-50 

E 

SIG(ll) 

atoms/cm 

atmospheric loss parameter Z 
for TIME(ll) 


51-60 

E 

SIG(12) 

atoms/ cm 

atmospheric loss parameter Z 
, for TIME (12) 

Card 6 

1-10 

E 

ENO 

# protons/ 
cm 3 

initial proton density 


11-20 

E 

TSUBO 

months 

initial time 

Card 7 

1-8 

F 

DTI 

months 

first cycle integration interval 


9-16 

F 

DT2 

months 

remaining cycle integration 
interval 


17-24 

F 

TEND 

months 

limit of integration 

Card 8 

1-3 

I 

IEST 

- 

initial energy level subscript* 


4-6 

I 

IEDEL 

- 

energy level subscript incre- 
ment* 


7-9 

I 

IEEND 

- 

final energy level subscript* 

Card 9 

1-5 

I 

IEPR 

- 

subscript of energy level to be 
printed on tape 5 


6-10 

I 

NPRINT 

- 

print on tape 5 after so many 
increments in time 


11-15 

I 

ICSUBO 


cycle corresponding to ENO 


and TSUBO 


see table of energy levels in INPUT. 
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GSFC FORM 541-1 (July - 60) 










H. Output 


Output from this program exists in several forms. First of all there is the 
output from the transient steady-state solution. In this solution there are two 
pages of printout for each case (each B-L line considered). The first page 
merely lists the input with explanatory headings. The second page lists each 
energy level and the cycle number at which steady-state was reached for that 
energy level. It also gives the maximum and minimum density and flux for that 
cycle with the times at which they appear within the cycle. Finally it lists the 
ratio of the maximum to the minimum flux. If, for any energy level, the absolute 
value of the minimum flux becomes greater than 50,000 protons/cm 2 sec. mev. 
then the program will ask that a smaller integration interval be used and go on 
to the next case. On the other hand, if the ratio of the maximum to the mini- 
mum flux becomes less than 1.09 for any energy level the program will neglect 
all higher energy levels and pass on to the next case since there will be no 
solar cycle variations in higher energy levels. This output occurs on tape 3. 

The time-averaged steady-state solution has the same output with the ex- 
ception that on the second page for a given B-L line $ , 1 , dE/dx and d/dE 
(dE/dx) are printed along with the flux and density for each energy level. 

Time histories for energy levels are printed on tape 5 if requested during 
the steady-state solution. Here time, flux and density are printed for each 
cycle at the intervals indicated by NPRINT. See figure 13 for an idea of the 
results to be expected from the transient steady-state solution. Figure 15 shows 
output from several time history runs. 
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Transient Steady -State Sample 
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TI^£ AVERAGE STEADY STATE SOLUTION FCK LATL * 15.945C , L = 1.2500 . B * 0.2240 
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Tape 5 Time History Sample 
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TRANSIcNl STEADY STATc SOLUTION FOR L AT D * 28.9050 » L * » H - D.22DB 
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Time Averaged Steady State Sample 
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THE HISTORY OF DENSITY ANO FLUX FUR CNfcRGY * 200.0 , lATO = 28 . 905 C , L = 1 . 6 UU 0 » 0 * 0.2208 
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I. Running Time 

Transient steady-state solutions where DTI = .5 and DT2 = 1.0 will cover 
eleven cycles per minute. See figure 14 for an idea of how many cycles are 
necessary to reach steady-state for various B's at an L of 1.25 earth radii. 

Runs requesting a time history print on tape 5 take about twice as long. 

Time averaged steady state solutions take about a quarter of a minute for 
each B-L line. 
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Vffl. STEADY-STATE CALCULATION 


A. Introduction 

This program evaluates the solar maximum and solar minimum steady 
state conditions of the conservation equation. That is, with dN p /dt set equal 
to zero, the flux and density are studied at solar minimum (time = 0.0 years, 
see Figure 1) and solar maximum (time = 4.0 years). The mean lifetimes of 
the protons are also calculated as well as dE/dx, d/dE(dE/dx), the three coef- 
ficients of the conservation equation C 0 , C x and C 2 and the source and loss 
terms. This is all printed together with the appropriate 2 for a given B, L, k 
and energy level. All 15 energy levels of the preceeding program are evaluated 
for each B and L. 


B. Equations 

The equations used in this program are listed below: 
SOURCE = XX = A a $/L 2 E B ° cos 4 k n 


A i N P /dE\ B d /dE \ B 

LOSS = YY = — + — (dx) + 


MEAN LIFETIME = TAU = N P /XX 

FLUX = FLUXP = N p (EXT)E B 2 C = N P /3C = N p v 


DENSITY = EN1 = 


A 0 $ 


L 2 E b o cos 4 A. 


,’ 4 (— ) 

°LE Bi \dX / 2 \dE \dX/ 


+ Z 


where C = speed of light 

$ = PREL = relative neutron source strength 
2 = SIG = atmospheric loss parameter 
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A Q , . . . , A = high or low energy conservation equation coefficients depend- 
ing upon whether E > 80 Mev. or E < 80 Mev respectively 

B 0 , . . . , B 2 = high or low energy conservation equation coefficients depend- 
ing upon whether E >80 Mev. or E <.80 Mev. respectively 


A. 0 = Mirror latitude 
N = proton number density 
v = neutron velocity 

fi = v/C 
EXT = /S/E® 2 


C . Mnemonics 


Quantity 

Description 

Units 

TIME(J) 

abscissa of time in increments of years for 
12 years 

years 

ELOSS(J) 

dE/dX corresponding to E(J) (see Figure 11) 

Mev. /cm. 

E(J) 

energy corresponding to ELOSS(J) 
(see Figure 11) 

Mev 

ALO 

A Q for E < 80 Mev. (see preceeding page) 

# protons 
cm. sec. Mev. 

AL1 

^ tl ft 11 

cm/ sec 

AL2 

^ U It ft 

cm. /Mev 2 sec. 

AH0 

A 0 for E > 80 Mev. (see preceeding page) 

# protons 
cm. sec. Mev. 

AH1 

^ 11 ft 11 

cm. sec. 

AH2 

^ H ff ft 

cm/Mev 2 sec. 
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Quantity 

Description 

Units 

BLO 

B 0 for E ^ 80 Mev. (see preceeding page) 

- 

BL1 

r> ♦» M Tf 


BL2 

g M 1 T H 

- 

BHO 

B Q for E > 80 Mev. (see preceeding page) 

- 

BH1 

O M f T It 

B i 

- 

BH2 

tj " It II 

B 2 


DELOSS(J) 

d(dE/dX)/dE corresponding to DE(J) 
(see Figure 12) 

cm.’ 1 

DE(J) 

energy corresponding to DELOSS(J) 
(see Figure 12) 

Mev. 

PREL(J) 

$ for HME(J) (see Figure 11 and 
page 167) 

- 

CONVR 

conversion factor - degrees to radians 

- 

ALAT 

mirror latitude \ 0 

degrees 

EL 

magnetic field line L 

earth radii 

B 

magnetic induction B 

gauss 

R(J) 

atmospheric scale factor R for TIME(J) 

- 

SIG(J) 

atmospheric loss parameter 2 for TIME(J) 

atoms/cm 

ALA TO 

ALAT in radians 

radians 

SVPR 

temporary storage of PREL(l) 

- 

SVSG 

temporary storage of SIG(l) 

atoms/cm 

t 

time 

years 
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Quantity 

Description 

Units 

M 

energy level subscript 

- 

ENER 

energy level E(M) 

Mev. 

ALOSS 

dE/dX for E(M) 

Mev/cm 

X(k) 

temporary storage of DE(k) 

Mev. 

Y(k) 

temporary storage of DELOSS(k) 

cm' 1 

EANS 

d(dE/dX)/dE for E(M) 

cm’ 1 

ALOSSA(J) 

(dE/dX) X R(J) for E(M) 

Mev/ cm 

DLOSSA(J) 

(d(dE/dX)/dE) X R(J) for E(M) 

cm’ 1 

AO 

temporary storage of ALO or AHO 

# protons 
cm. sec. Mev. 

Al 

temporary storage of AL1 or AH1 

cm. /sec. 

A2 

temporary storage of AL2 or AH2 

cm/Mev 2 sec. 

BO 

temporary storage of BLO or BHO 

- 

B1 

temporary storage of BL1 or BH1 

- 

B2 

temporary storage of BL2 or BH2 

- 

EXT 

(see page 168) 

- 

CO 

conservation equation coefficient c Q where 
C Q = A q /L 2 E b ° cos 4 ^ 0 

# protons 
cm 2 sec. Mev, 

Cl 

conservation equation coefficient C where 

C = A /E B i 
1 

cm/ sec. Mev. 

C2 

conservation equation coefficient C 2 where 

C = A E b 2 
2 2 

cm/sec. Mev. 
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Quantity 

Description 

Units 

KCK 

counter to determine if solar maximum or 
solar minimum is being evaluated 

- 

T1 

particular time under consideration 

years 

EN1 

proton density 

# protons 
Mev. cm 

FLUXP 

proton flux 

# protons 
cm 2 sec. 

TAU 

mean proton lifetime 

sec 

XX 

conservation equation source term 
(see page 167) 

# protons 
cm 2 sec. 

YY 

conservation equation loss term 
(see page 167) 

# protons 
cnPsec. 
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D. Flow Charts 


1. Main Program 


START 

















2. Subroutine Table 


START 
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EXTERNAL rCRMILA NUMBERS With CCH Ku SPCNC I NG INTERNAL FURKULA NUMBERS AN L OCTAL LOCATIONS 
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50 IF(LL) 41,41,4? 

41 YANS=ANS 
GO TO 2 3 

42 YAN$ = tXPF ( ANS ) 
23 RETURN 
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F . Input 


Input to this program occurs on logical tape 2 and consists of B, L, k, R(l), 
..., R(12) 2 (1 ),...,£ (12) for each case under consideration. As with the con- 
servation equation calculation the B-L- 2card is followed by two R cards which 
are followed in turn by two 2 cards. The control cards are left out of this pro- 
gram since there is only one method of computation available here. 


1. Input Card Description 



Columns 

Mode 

Quantity 

Units 

Description 

B-L-A. 

1-8 

F 

ALAT 

degrees 

latitude 

Card 

9-16 

F 

EL 

earth 

radii 

magnetic field line 


17-24 

F 

B 

gauss 

magnetic induction 

First 

R Card 

1-10 

E 

R(i) 

- 

scale factor R for TIME(l) 


11-20 

E 

R(2) 

- 

scale factor R for TIME(2) 


21-30 

E 

R(3) 

- 

scale factor R for TIME(3) 


31-40 

E 

R<4) 


scale factor R for TIME(4) 


41-50 

E 

R(5) 

- 

scale factor R for TIME(5) 


51-60 

E 

R(6) 

- 

scale factor R for TIME(6) 

Second 

R Card 

1-10 

E 

R(7) 


scale factor R for TIME (7) 


11-20 

E 

R(8) 

- 

scale factor R for TIME (8) 


21-30 

E 

R(9) 

- 

scale factor R for TIME(9) 


31-40 

E 

R(10) 


scale factor R for TIME(10) 


41-50 

E 

R(ll) 


scale factor R for TIME (11) 


51-60 

E 

R(12) 

- 

scale factor R for TIME(12) 

First 

2 Card 

1-10 

E 

SIG(l) 

atoms /cm 

atmospheric loss parameter 2 
for TIME(l) 


11-20 

E 

SIG(2) 

atoms/cm 

atmospheric loss parameter 2 
for TIME(2) 


21-30 

E 

SIG(3) 

atoms/ cm 

atmospheric loss parameter 2 
for TIME (3) 


31-40 

E 

SIG(4) 

atoms/ cm 

atmospheric loss parameter 2 
for TIME (4) 
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Columns 

Mode 

Quantity 

Units 

Description 

41-50 

E 

SIG(5) 

atoms/cm 

atmospheric loss parameter £ 
for TIME(5) 

51-60 

E 

SIG(6) 

atoms/ cm 

atmospheric loss parameter Z 
for TIME (6) 

Second 

I Card 1-10 

E 

SIG(7 ) 

atoms/ cm 

atmospheric loss parameter Z 
for TIME (7) 

11-20 

E 

SIG(8) 

atoms/ cm 

atmospheric loss parameter Z 
for TIME (8) 

21-30 

E 

SIG(9) 

atoms/cm 

atmospheric loss parameter Z 
for TIME(9) 

31-40 

E 

SIG(10) 

atoms /cm 

atmospheric loss parameter Z 
for TIME (10) 

41-50 

E 

SIG(ll) 

atoms/cm 

atmospheric loss parameter Z 
for TIME (11) 

51-60 

E 

SIG(12) 

atoms/ cm 

atmospheric loss parameter Z 


for TIME (12) 
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Sample GENERAL PURPOSE DATA SHEET 
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CSFC FORM 541-1 (July - 60) 



G. Output 


Output for this program occurs on tape 3. Time, density, flux, mean life- 
time, dE/dX, d(dE/dX)/dE, 2, source and loss are all printed for each energy 
level at a given B, L and latitude. The headings are listed below with explana- 
tions and dimensions: 

TIME — time — years 

N(E) — density — # protons/cm 3 Mev. 

FLUX - flux — # protons/cm 2 sec. Mev. 

MEAN LIFETIME - proton mean lifetime — sec. 

CO - C Q (see MNEMONICS) — # protons/cm 3 sec. Mev. 

Cl - C t (see MNEMONICS) - cm/sec. Mev. 

C2 - C 2 (see MNEMONICS) - cm/sec. Mev. 

ALOSSA - (dE/dX) ATMOS - Mev./cm. 

DLOSSA - (d(dE/dX)/dE) ATMOS - cm' 1 

SIG - I - atoms/cm 

SOURCE — conservation equation source term — # protons/cm 3 Mev. 
LOSS — conservation equation loss term — # protons/cm 3 Mev. 

ENERGY - energy - Mev. 

LA TO - latitude - degrees 
L — magnetic field line — earth radii 
B - magnetic induction - gauss 
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L ATO* 


8 = 0.2331 


12.954 L - 1.18b 
ENERGY- lO.QOMEV 


TIME 

N( t) 


FLUX 

MEAN LIFETIME 

0. 

0.71 726-13 

0.00031 

0.306936-04 

CO* 0 .72115 E— 15 

Cl* 0.10386E 09 

C2* 

0.21759E LO 


.1 >34* 0.47404E-10 

D L 0 S S A = 0. 35124 E ”11 

SIG* 

0. I3O80E-14 


nIF> = 0 . 7 1 7 2 0 E - i 3 

SOURCE = 0.90144E-15 

LOSS* 

0. 12 569 C — 01 


TIME 

N { E ) 


FLUX 

MEAN LIFETIME 

o 

6 

• 

0. 70836-13 


0.00031 

0 . 3 7 8 9 3 F “ 0 4 

CO* 0.721 1 5E- 1 5 

Cl* 0.103866 09 

C 2 * 

0.21 / 3 9 E 10 


AL OSS A- 0.384006-10 

DLOSSA* 0. 28452 E“ 11 

SIG* 

0* 1 0 6 1 0 c - 1 4 


N(£) = 0.70831E-13 

SOURCE* 0.72115E-L3 

LOSS* 

0. 101316-01 


fc: N £ R GY = 2‘j.OOMbV 




1 IME 

NIC) 


FLUX 

MEAN LIFETIMt 

0. 

0.2746E-13 


0.00013 

0.95765E-04 

CO* 0.72376E-16 

Cl* U.64318E 08 

C2* 

0 . 3 3 5 8 6 C 10 


ALQSSA- 0.22595G-10 

OlUSSA* 0.76320E-I2 

SIG* 

0. 130806-14 


N(E) = 0.22457E-13 

SOURCE = 0.90470E-16 

LOSS* 

0. 40286C-02 


TIME 

NIE) 


FLUX 

Me. AN L I FET I ML 

4.000 

0.2213E-13 


0.00015 

0. 1 1822E-03 

Co = 0.72376E-I6 

C 1= 0.64313C 03 

C 2 * 

0.3368GC 10 


ALOSSA* 0 • 133036*10 

DLOSSA* 0.61B23E-12 

SIG* 

0. 10ol0fc-i4 


N ( F ) = 0.22178E— 13 

SOURCE* 0 . 72 3 76 c- 1 6 

LOSS* 

0 • 32634c — 02 


ENERGY* sO.OOMEV 




T IME 

NIE) 


FLUX 

MEAN LIFETIME 

0. 

0.998 IE ”14 


0.00OU9 

0 • 2 4 2 2 7 E - 0 3 

CO* 0.12715E-16 

Cl* 0.44 760 E 03 

C 2* 

0 • 4 6 b 8 6 1 10 


A L 0 S S A - 0. 12956G-10 

OLUSSA* 0*2146 3 E”12 

SIG* 

0. 13U80E-14 


N { E ) = 0.99806E-14 

SOURCE- 0. 13894E-16 

LOSS* 

0. 159256-02 


TIME 

N ( F ) 


FLUX 

MEAN LIFETIMt 

4.000 

0 . 9 8 3 7 C ” 1 4 


0 • 0 0 0 U 9 

0. 2990b E-03 

CO* 0. 12715E-16 

Cl* 0.44 760E 08 

C2* 

0.46886E 10 


ALOSS A * 0* 10 4 95E- 1 0 

OLUSSA* 0.17387E-12 

SIG* 

0. 10ol0b-14 


N { F ) = 0.983676-14 

SOURCE* 0.12715E-16 

LOSS* 

0. 12900E-02 


ENERGY* 

ri.OOMEV 




TIKE 

N { F ) 


FLUX 

MEaN LIFETIME 

0. 

0.63o8c-14 


0.00007 

0 . 42 73 3 E— 03 

CO- 0.45973E-17 

Cl* 0.36207E 08 

02 = 

0.568916 10 


ALOSSA* 0 . 9 4 4 3 I F - 1 l 

DLOSSA* 0 * 9 7 2 1 4 E - 1 3 

SIG* 

0. 130806-14 


N(F) = 0.636816-14 

SOURCE* 0 . 5 7 4 6 6 E “ l 7 

LOSS* 

0. 902416-03 


TIME 

M ( 6 ) 


FLUX 

MlaiN LIFETIME 

4.000 

0.62896-14 


0.00007 

0.527776-03 

CG = 0.459736-1 7 

Cl* 0.36^07F 03 

C 2 = 

O.Sbb^iL 10 


ALOSSA* 0 . 7 6 4 9 3 E - i 1 

DLOSSA* 0.78/46E-13 

SIG* 

0. 106 l OE- 1 4 


M(G)= 0 • 628 90l- i 4 

SOURCE* 0. 459736-17 

LOSS* 

0. 73 l 0 0£ -03 


ENERGY* 1G0.00M6V 




TIME 

N ( E ) 


FLUX 

M £ A N LIFETIME 

0. 

0.48 9 3E- L 4 


0.0000b 

0 . 69?44E”0 3 

CO* U. 2 1809c- I 7 

Cl* 0.22^096 Ob 

C 2 * 

0.65h75E 10 


ALOSSA* 0. 76002b- 11 

DLOSSA* 0.5 76596— I 3 

SIG* 

0. 13080E-14 


N ( C ) = U • A8928C - i A 

SUUKCc* 0. 27261 E-i 7 

LOSS* 

0.567166-03 


TIME 

N 1 t ) 


FLUX 

MEAN LIFETIMt: 

A. 00 0 

0.4832E-14 


0.0000b 

0.85460E-03 

CO = 0.21 a09E- 1 7 

Cl* 0.22 :>09E 08 

C2* 

0. 63475 t 10 


ALOSS A= 0. 6 1 56SC- 1 I 

DLOSSA* 0.467076-13 

SIG* 

0. 10olGc-14 


N ( C ) = 0. 48 A20E- 1 A 

SOURCE* 0.21 5096-17 

LOSS* 

^ 0.431 34E-03 


E 'IE <GY = 123.00MEV 




TIME 

NIE) 


FLUX 

MEAN LIFETIMt: 

0. 

0.38 12c- 14 


0.00003 

0.952896-03 

CO = 0.12348E-17 

Cl* 0.194446 U8 

C2 = 

0 • 7 G o 9 9 E 10 


ALQSSA* 0.6459 9 E - 1 I 

DLOSSA* 0.381936-13 

SIG* 

0. 130806-14 


N(6) = 0 • 381236-14 

SOURCE * U. 15435F-17 

LOSS* 

0* 40488E ”03 


\ IMF 

NIE) 


FLUX 

MEAN LIFeTIMt 

4.000 

0. 376 5 G — 14 


0.00005 

0. 11763E-G2 


186 


C0= U.12348E-17 

Cl= 0.19444E 08 

C2 = 

0. 705996 l o 


ALOSSA= 0. 523286-11 

OLDS S A= 0 . 309 3 HE — 13 

S I G- 

0. lUulOE-14 


N ( E ) = U. 37650E-1 4 

SOURC£= 0.12348E-I7 

LUSS = 

0 • 32 798C-0 3 


ENERGY 3 I50.00KEV 




TIME 

N{£) 


FLUX 

MEAN LIFETIME 

0. 

0.31 13E-14 


0.000U5 

0.123916-02 

CO 3 0 . 7 7 5 84E- 1 8 

C 1 = 0.172526 08 

i/2 = 

0.732756 10 


A $SA = 9.56o31E-ll 

OLOSS A- 0.2 702 9£- 1 3 

S I G = 

0. 13080E-14 


W<C) = U.31148E-14 

SGURCE= 0. 96 98 OE- 18 

LOSS = 

0. 31135E-03 


TIME 

N ( F ) 


FLUX 

MEAN LIFETIME 

A. 000 

0.3076E-14 


0.00005 

0. 152966-02 

C0= U.77584E-18 

Cl- 0.17232E 08 

C2~ 

0. 732756 10 


ALOSS A = 0. 46035E-1 1 

DLOSS A= 0.218956-13 

S 1 G = 

0 . 106 10E— 14 


N ( E ) = 0.30/61E-14 

SOURCE^ 0.77bo4E-18 

LOSS = 

0.23222E-03 


ENEKGY= 200 . OOMb V 




TIME 

NIE ) 


FLUX 

MEAN LIFETIMt 

0. 

0.2233E-14 


0.00004 

0. 18 701E-02 

CO- 0.3726BE-18 

Cl= 0.142B5E 03 

C2 = 

0.83106 6 10 


ALOSSA 3 0.46921E-11 

L)LOSSA= 0.154 51c— 13 

S 1 G = 

0. 13080L-14 


NIC) 3 0.22579E-14 

SOURCE 3 0 . 46 j32 E-18 

LOoS = 

0 * 2Ut>30E — 03 


TIME 

MIC) 


FLUX 

MEAN LIFETIMt 

4.000 

0.2230E-14 


0.00004 

0.230846-02 

CO 3 U . 37265E- 1 o 

Cl= 0.14285E 08 

C2 = 

0.83106E 10 


A L 0 S S A = 0.3800BE-11 

OLUSSA= 0. 12516E— 13 

S I G = 

C. 106 IOC- 14 


N(E)= 0 . 22298E— 14 

SOURCE 3 0. 37265E-1B 

LOSS = 

0. 16/13E-03 


F N E K G Y = 230 . DOME V 




r imf 

NIC) 


FLUX 

MEAN LIFETIME 

0. 

0.1733E-14 


0. 00003 

0.256416-02 

C0= U.21100£-iB 

Cl 3 0.12 a40E 08 

C2 = 

0.89736E 10 


ALOSS A = C.40869E-11 

DLOS S A= 0 . 98 392E-14 

S I G = 

0.130806-14 


N I E ) = 0.17529E-14 

SOURCE 3 0. 26375 E- 18 

LDSS = 

0. 13046E-03 


TIME 

N ( E ) 


FLUX 

MEAN LIFETIME 

4,000 

0. 1731E-14 


0.00003 

0.31651 £—02 

C0= 0.21100E-18 

Cl= 0.12340E 08 

C2 = 

0.89/366 10 


AL OSS A = 0.33106E-11 

OLOSSA 3 C.79702E-14 

S I G= 

0. 106l0t-14 


N { 6 ) = 0* 17310E-14 

SOURCE 3 0. 21 IOOE-I8 

LOSS = 

0.121 89E— 03 


ENERGY 3 300.00MEV 




TIME 

N{E) 


FLUX 

MEAN LIFETIME 

0. 

0. 1404E-14 


0 • OOOu 3 

0.326806-02 

C0= 0. 13257E-18 

Cl= 0.10949E 08 

C2 = 

0.933456 10 


ALOSS A- 0.36805E-I1 

DLOSSA= 0.683036-14 

S1G- 

0.130806-14 


N { t ) - 0* 140376- 14 

SOURCF* 0.163716—18 

LOSS = 

0.116056-03 


TIME 

N ( E ) 


FLUX 

MEAN LIFETIMt 

4.000 

0.1386E-14 


0.000U3 

0.403386-02 

C0= 0.13257E-16 

C 1 = 0.10349E 08 

C2 = 

0.953456 10 


ALOSS A = 0.29813E-11 

DLOSSA= 0.55329E-14 

S 1 G = 

0. 106106-14 


N ( E ) * 0. 133616-14 

SOURCE = 0.13257E-18 

LOSS = 

0 • 93o43£ -04 


E\EKGY= 350.00MEV 




TIME 

N (E) 


FLUX 

MEAN LIFETIME 

0. 

0.1142E-14 


0.00002 

0-393 76E— 02 

C0= 0.894956-19 

C 1= 0.98958E 07 

C2 = 

0.10U75E 11 


AL OSS A= 0.33696E-11 

DLOS S A = 0.50877E-14 

S I G = 

0.13080E— 14 


N ( E ) = 0.1 14186—14 

S0URC6= 0.11 187E— 18 

LOSS = 

0. 979 78E“04 


TIME 

N ( E ) 


FLUX 

MEAN LIFETIME 

4.000 

0 . 1 l 27£- 14 


0-00002 

0 • 4860 1 E— 02 

C0= 0. 89495E- 1 9 

C 1 = 0.98958E 07 

C2 = 

0.10075E 11 


ALOSS A= 0.27458E-11 

DLOSS A= 0.41212E-14 

S IG = 

0 • 106 10E — 14 


N ( E ) = 0.112 f4E— 14 

SUURC£= 0.89495E-19 

LOSS = 

0. 79331E-04 


ENERGY^ 400.00MEV 




TIME 

Hit ) 


FLUX 

MEAN LIFETIME 

0 . 

0.9561E-15 


0.00002 

0 • 46341 E— 02 

C0= 0. 63676E-19 

Cl= 0.906596 07 

C2- 

0.103486 11 


ALOSS A= 0.31 723E- 1 1 

DLOSS A= 0.38580E-14 

S I G= 

0- 13080E-14 


N(E)= 0.95607E-13 

SOURCE= 0.79596E-19 

LUSS= 

0 • 8325 3E-04 


TIME 

Hi E) 


FLUX 

MEAN LIFETIME 

4.000 

0.9440E-15 


0.00002 

0.57195E-02 
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CO* 0# 636 76E- 19 

C 1 = 0. 90659L 07 

C2 = 

0 . i 0 D 4 8 L 11 



ALflSS A= 0. 25697F-1 1 

OLOSSA* 0 . 3 1 2 5 2 E - 1 4 

S l G = 

G * 10ol0t-14 



N(E) = 0.94A00E-15 

SOURCE* 0.636 76 £— 19 

LOS S = 

0.6 74 54b -04 



energy* soo.oo’iev 





T I Mb 

Nit) 


FLUX 

MEAN 

L IFETIMt 

0. 

0 . 7_>9 f> 3£ - 1 5 


0.00U02 

0. 595 1 8E-C2 

CO = Q.*fco r >4E-19 

Cl* 0.7831 3E U 7 

C2 = 

0.11390b 11 



SS^ = D.26715E-11 

OLOSSA* 0 . 24088c - l A 

S l G = 

0. 13U60c-l4 



N(c) = 0.69626E-19 

SJURCE* 0. A 506 BE- 19 

LOSS = 

0.64d2iE-04 



TIME 

NIL ) 


FLUX 

Mi E A N 

LiFETIMh 

4.000 

0.6864E-10 


0.00002 

0.73451b -02 

C0 = 0.36054F-19 

Ct = 0. 78313E 07 

C2 = 

0.11390E 11 



AL OS S A = 0.?3260F-ll 

OLOSS A= 0. 19512c- U 

S 1 G = 

0. 10610b- 14 



M ( E ) = 0 .68 642 1 - 1 6 

SOURCE* 0 . 360545-19 

LOSS* 

0. 5232 6 E- 04 



C'lcKGY* 6 0 0 . 0 OM £ V 





T IMt 

M { !r > 


FLUX 

MEAN 

LIFETIME 

0. 

0 . 52486— 1 5 


0.00001 

0.71 4 99 E— 02 

C0= 0.22653E-I9 

Cl* 0.69485E 07 

C2 = 

0.12127E II 



ALOSS A= i; . 26 76 IE- 1 1 

OLOSSA* 0. 16081E-1A 

S I G = 

0. 13080E-14 



N(E) = 0 . 6 2 4 7 7 L - 1 5 

SOURCE* 0.26316F-19 

LOS S = 

0. 539596-04 



TIME 

NIE ) 


FLUX 

MEAN 

LIFETIME 

A. 000 

0.51815-15 


O.OOOul 

0.882 30E-02 

C0= 0.22653E-19 

Cl* 0.694856 07 

CP- 

0.1 212 7E 11 



AL0SS4* <J.2lo78f-ll 

ULQSS A= 0.13026E-1A 

Sib* 

0. 1 Obi Ob- 14 



M(E)= 0.51 <i 0 •"> E — 1 5 

SOURCE* U. 22653 E— 19 

LOS S = 

0.43727 l -04 



ENERGY* 7U0.00MEV 





T I ME 

MID 


FLUX 

MEAN 

LIFETIME 

0. 

0.4033H-15 


0.00001 

0 . 6 1 39 96-G2 

C0= 0. 16292F-19 

Cl* 0.62 602 E 07 

02- 

0. 12788b 11 



A L 0 S S A = 0 * 254146— i 1 

OLOSS A= 0 . 1 1 50 3E- 1 A 

SI G= 

0. 1 3 0 8 0 b - 1 4 



iM ( E ) — 0.40 3311-19 

SOURCE* 0. 19115 E- 19 

LOSS* 

0.47397E-04 



TIME 

0 ( - ) 


FLUX 

MEAN 

LIFETIME 

4.000 

0.3981E-15 


0.00001 

0. lOOAAfc— G 1 

C0= G.15292L-19 

Cl* 0.628C2E 07 

C 2 = 

0.12/88E 11 



ALOSS A= 0 • ?0 5 86E- i 1 

OLOSSA* 0.93180E-15 

S I G = 

0. 106106-14 



N ( E ) - 0. 3 9 oil £-13 

SOURCE* 0.152926-19 

LQSS = 

0. 3o4 1 2E-04 
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H. Running Time 

This program will do eleven cases in a minute and a half. 
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REFERENCE TIME, TO 
TIME 

Figure 1 -A time history of the 10.7 cm solar flux according to the measurements of the National Research Council 
of Canada for the recent past. The heavy dotted line indicates the approximate yearly average. 
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Figure 2-A mapping of the polar coordinates R and K onto the B-L plane 
where R = (h +6378.2) 
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Figure 4-B-L contours for the southern hemisphere at an 
L of 1 .25 earth radi i 
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WHERE: 

d$ - Element of Arc along the particle's helical trajectory 
dl - Element of Arc along the field line 

Figure 5-Schematic of a trapped particle’s north-south motion 
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Figure 6-The weighing factor, A(A), versus latitude for various mirror 
latitudes, A , where A Q are the asymptotes of each curve 
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Figure 13-A comparison of the steady-state and transient proton flux energy spectrums 
for L = 1.25, B» = .199 at solar minimum and solar maximum 


205 






Figure 15-A time history of proton flux as a function of energy at L - 1,25, 
B = ,209 for the first and tenth solar cycles 
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